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I INTRODUCTION 


A. CFAR SIGNAL PROCESSING 


Signals processed by radars are inherently corrupted by environmental clutter and 
thermal noise from receiver components. The echo signal returned from any potential target 
will thus be accompanied by noise from these sources. The simplest method of detecting 
the target return signal is to set a fixed power amplitude threshold level which is known to 
be greater than the mean noise level and most of the noise peaks, but less than the expected 
magnitude of the target return signal. Thus, any signal which is measured to be greater than 
the threshold is declared by the radar to be a target detection, as shown in Figure 1. Such a 
detection may include a noise peak, and thus a false detection. This is known as a false 
alarm. By measuring the average number of false alarms occurring in a specified period of 
time, a false alarm rate can be established. The ratio of the average number of false alarms 
to the total number of noise samples is known as the probability of false alarm (Pra). By 
the same token, the ratio of the average number of target signal detections to the total 
number of target signals sampled is known as the probability of detection (Pp). This 
implies that some target signals may not be detected, and thus the probability of a missed 
target detection can be computed. In general, the detection performance of a radar is stated 


in terms of the Pg, and the Pp. 


Fixed threshold detection is inadequate for radar systems that have to detect the 
presence of targets within background environments that are more complex and less known 
than thermal noise alone. This type of detection produces large changes in false alarm rates, 
a very undesirable characteristic for radar systems intended for high target detection 
sensitivity and reliability. Also, the great number of detections that occur saturate and 


disrupt the tracking computer associated with the radar system. Thus detector design 


criteria for modern military radars is to maintain a low and constant Ppa as well as a good 


Pp for a given signal-to-noise power ratio (SNR). Automatic detection systems which 
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Figure 1: Detection Threshold 


maintain a constant Ppa are known as constant false alarm rate (CFAR) systems. In these 


systems, the constant false alarm rate 1s maintained by adaptively varying the detection 
threshold level in accordance with estimates of the noise power (targets not present). Such 
detectors can easily cope with fluctuating noise levels, as the estimates are taken in real 
time. CFAR detection is usually performed in range. That is, the noise estimates are taken 
from a predetermined number of range bins (also known as cells) which are adjacent to the 
range bin on which detection is being performed (known as the test cell). This is shown in 


Figure 2. These adjacent cells are called reference cells. The analog-to-digital converter 


samples taken over the time period occurring between two consecutive transmitted pulses 


are the range bins. That is, each range bin represents an increment in range from the radar. 


transmitted radar pulses 


test cell 


reference cells 


Figure 2: Unambiguous Radar Range with a Sample CFAR Window 


The maximum possible range for the time period between the two pulses is known as the 
unambiguous range. The reference cells and test cell form a window which slides through 
a specified set of range bins, evaluating each range bin for a target. 

Most adaptive threshold detectors assume that the samples in the reference range cells 
are independent and identically distributed. Furthermore, it is usually assumed that the time 
saraples are independent. In essence, the detectors test whether the test cell power is 
sufficiently larger than the noise power in the reference cells. It is assumed that the 
probability density function of the noise is known except for a few unknown parameters. 
The surrounding reference cells are then used to estimate the unknown parameters, and a 
threshold based on the estimated parameters is obtained. Since the determined parameters 


are estimates, the resulting threshold has some error and must therefore be slightly larger 


than the threshold that one would use if the parameters were known exactly a priori. This 
causes a loss in target detection sensitivity, called a CFAR loss, and thus results in a 
degradation of the Pp. 

In order to analyze the Pp and the Pg, performances theoretically, the probability 
density functions of the radar range bins have to be derived. In general, given a density 
function for a non-target cell (noise only) and one for a target cell (signal with noise), the 
Ред and the Pp are obtained by integrating over these functions for a given detection 


threshold T, as shown in Figure 3. Obtaining these functions can be very tedious. The 
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Figure 3: Pp and Pra 


resulting analytical expressions representing these density functions are not always 


practical to manipulate and program. Thus, attempts are often made to model the equations 


with simple closed form expressions to facilitate the analysis and implementation of the 


desired CFAR processing scheme. 


B. PRINCIPLE CONTRIBUTIONS 

The main contribution of this research is the rigorous treatment of the effects of two 
envelope detection approximations on the statistical characteristics of the detected signals 
within a greatest-of CFAR (GO CFAR) signal processor. The first approximation is given 
by x, = aMAX {|/|,|Q|} + bMIN {|1, |Q] } , where I and Q are, respectively, the in-phase 


and quadrature samples of the input signal. The values a and b are simple multiplying 


constants. A closed form solution for the probability density function of this detector is 


shown to be untrackable analytically. An alternative method to obtain a closed form 
expression for the density function is developed. This approach curve-fits the results of a 
numerical convolution. These closed form expressions enable a quick implementation of 
the function within system simulation programs. Reference tables provide the necessary 
data for this method. 

The second approximation investigated is given by X, = a|/| +b|Q|, a simplified 
form of x,. An analytical closed form solution is presented, as well as a curve-fit solution 
similar to that obtained for x,. This allows a comparison of the two envelope detection 


approximations in terms of their coefficients. 


Finally, a comparison between the true envelope detector GO CFAR performance and 
the performance using x, and £, is made in terms of numerically derived probability 
density functions and the Pp. Results show that there is a significant difference in 


processor performance for the case of / and Q being zero mean white Gaussian noise with 


a non-fluctuating target. 


C. THESIS OUTLINE 


The objective of the research documented here is to investigate a method of 


computing the Pp analytically for a GO CFAR processor that uses the envelope detection 


approximations x, = aMax {|/\,|Q\} +bMin{|\,|Q|} and x, = al/| + b|Q|. The main 
hurdle in doing this is to determine a closed form expression for the probability density 
function of a detected radar range cell with a target signal present in it. 

Section II describes the operation of a GO CFAR in a simple manner. The 
functionality of it is described by stepping through all its operations sequentially, providing 
a brief discussion at each point. It then provides a short discussion of envelope detection 
approximation. Finally, the requirements for the analytical determination of the Pp for the 
GO CFAR detector is outlined. Section III provides the rigorous analytical development of 
the closed form solution for the probability density function of a target test cell using the 


envelope detection approximation x,. It is a step-by-step walk-through of the derivation, 


presented first as a general result and then applied to the particular requirements of this 
research. The difficulties encountered and the assumptions made for every step of the 
derivation are described. This section also provides a closed form expression derived from 
curve-fitting the numerical convolution. Section IV presents the closed form solution for 


the probability density function of a target test cell using the envelope detection 
approximation x,. A closed form solution using curve-fitting techniques is also described. 
Section V presents a comparison of the two methods of envelope detection approximation 
discussed in Sections III and IV in terms of the Pp and the Pra. Section VI provides 
concluding remarks on the success of this research. It also discusses the possibilities for 
continued research in this field, as well as different aspects of the problem which could be 
further studied. 


D. REFERENCE OVERVIEW 


The studies by P.E. Pace and L.L. Taylor [Ref. 1] concentrate on obtaining Ше Ред 
analytically for an envelope detection. approximation GO CFAR processor. The 


approximation used is X, — al|/| -- b|Ql. They derive a closed form solution of the 


probability density function of a radar range cell with noise only in it. Being able to 


integrate this function in closed form, they are able to produce Pra curves numerically. The 
Pra curves are derived as a function of a threshold multiplier and the number of reference 
cells used on each side of the test cell. They also derive, through curve-fitting, closed form 
expressions for the Pg, and compare its performance to the numerical solution. 


The technical memorandum by D.J. Wilson [Ref. 2] derives the Pp for a test cell with 


an envelope detection approximation £, with a = 1 and b = 1. He provides a useful insight 
into the analysis of this type of approximation and the subsequent derivation of the Pp. 

Another technical memorandum, this time by D.G. Loberger [Ref. 3], compares the 
performance of a cell-averaging CFAR processor (similar in concept to GO CFAR 
processing) for two different configurations. He analyzes the situation in which the test 
range cell is not used in the noise averaging of the reference cells. He also examines the 
situation in which the test cell is included in the averaging. The important result of this 
paper is to relate the threshold multipliers of each situation considered with a simple 
conversion formula. Using this result as a constraint guarantees the same Pra and Pp 
performances for either of the envelope detection approximations. 

A study by G.P. Laulusa [Ref. 4] evaluates the Ppa and Pp performances of a GO 
CFAR processor which includes the test cell along with the reference cells in the estimation 
of the noise level. He uses the envelope detection approximation £, with unity multiplying 
coefficients and a normalized version of the probability density function of the test cell. 

The paper by V.G. Hansen and J.H. Sawyers [Ref. 5] determines the CFAR loss of a 
GO CFAR scheme as compared to a simple cell-averaging CFAR. In the process of doing 
this they derive an integral expression for the Pp in terms of the probability density function 
of a greatest-of output block. This, in effect, defines the main building block for any Pp and 
Pfa computations on a GO CFAR processor, for any type of signal detector at the input of 


the processor and for any variation of the GO CFAR processing scheme. 


П. GO CFAR OVERVIEW 


A. IMPLEMENTATION 
A typical GO CFAR processor is shown in Figure 4. It uses an envelope detector, a 
lead window of reference cells, a lag window of reference cells, a test cell, two adders, a 


maximum value detector, two multipliers and a binary comparator. The true envelope 


detector is x — J + Q2, where / and Q represent the in-phase and quadrature samples, 
respectively, of the input signal. . The / and Q samples are obtained from an analog-to- 
digital converter (ADC), which produces the digital representations of the radar return 
signal. The ADC receives its analog / and Q signals from the synchronous detection and 
down conversion circuits. 

The envelope detected signal samples (from the / and Q samples) are not actually 
written into a GO CFAR processor per say, but rather are stored in consecutive range cells 
(or bins). This is not depicted in Figure 4. The GO CFAR processor is a sliding window 
which reads in the signal samples of the necessary range cells for processing. The sliding 
window enables the GO CFAR processor to test each cell for potential targets. It is divided 
into lead and lag segments, each of length n cells. Selection of n depends on the desired 
quality of the noise estimate, and usually lies between 4 and 32. 

The signal samples of the lead and lag windows are summed and the two results 
compared. The maximum sum is selected to calculate the estimate of the noise level. This 
value is then multiplied by a factor of 1/n to provide the final time-average noise estimate 
for the signal in the test cell. The estimate is then multiplied by T, a threshold multiplier, to 
establish a final detection threshold V,. This multiplier 1s preset to a value which will yield 


a desired, constant Pra. Finally, the detection threshold is compared to the signal present 


in the test cell. If the signal in the test cell is larger than the detection threshold, a target 
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Figure 4: GO CFAR Processor 


detection is declared. It can be concluded that GO CFAR processing is fairly simple in 


principle. However, it is much more complex to analyze theoretically. 


B. ENVELOPE DETECTION APPROXIMATION 
The envelope detector can be implemented with the envelope detection 


approximation mentioned in section I of the form 


x, = aMax (വ, 01) 222 (0, 01) . (1) 


Another similar approximation is 


7, с ад РО. (2) 


The constants a and b of the envelope detection approximation are predetermined so as to 
yield a desired performance level for the approximation given in (1). Envelope detection 
approximation is performed mainly for two reasons. First, much less hardware is required 
for implemention (less expensive). The squaring operation and the square-root operation 
in the true envelope detector are also difficult to design digitally so as to yield answers with 
minimal numerical rounding errors. 

Secondly, envelope detection approximations perform faster in real time compared to 
the true envelope detector. This is a crucial requirement, since the sequential GO CFAR 
processing on all the range cells must be completed prior to the next received pulse. For this 
reason, GO CFAR processors are frequently found in low pulse repetition frequency 
(LPRF) radar modes of operation. 


The Max and Min operators in the envelope detection approximation x, have the 


10 


effect of uniformly distributing each complex input sample, of the form Š = / + jQ, over 
the range of angles 9 = 0... (п/4) radians. This dependence on 6 can be expressed by 


rewriting the Max and Min terms in polar coordinates as 


፲ = 2, (2- 050 3-7 - 520) | (3) 


The relative approximation error can be defined as 


X, — X 


Ле 


е (9) = = 1- (а: со509 +В: 5110) ` (4) 





Then, the performance criteria in terms of the average error and mean-square error can be 


expressed as, respectively, 


2 
0 
and 
Цоо 
е? = - | е? (0) 40 . (6) 


A number of envelope detection approximations are shown in Table 1. [Ref. 1]. 

The coefficients of Table 1 are considered in this research. This allows the findings 
from this thesis to be correlated to the results of previous research in this area [Ref. 1]. They 
provide a good variety of performances, in that the errors can be observed to vary over a 


wide range. The sixth approximation is such that it yields a zero average error. 
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TABLE 1: ENVELOPE DETECTION APPROXIMATIONS 


approxi- 
mation 
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C. ANALYTICAL DETERMINATION OF PROBABILITY OF 
DETECTION 


It has been assumed for this paper that the / and Q input samples to the GO CFAR 
processor are independent and identically distributed. It is also assumed that the noise 
components of the samples are normally distributed zero mean and unity variance 
(~N(0,1)). Once the probability density function of the target test cell has been obtained, 


the probability of detection Pp can be determined using 


Pp(SNR) = 2 [ Py, n (2) ІШЕ d | | Pc (x, SNR) ' dz . (7) 
0 0 Tz 


n 


The probability of false alarm Pra is given by [Ref.1] 


Pra = 2[Py,n(2) Шо g | [р.б « doo (8) 
0 0 Tz 


The term p, ,, represents the density function for the n reference cell summation and р, 
represents the density function of the test cell. From Figure 4, መ. is Obtained by 
performing an n-fold convolution of the density function of a single reference cell. That is, 
each reference cell is assumed to have the same noise distribution and to be independent of 
each other. The distribution of the term p, is the only difference between the Pp (SNR = 0, 
target present) and the Pra (SNR = O, no target, noise only). In the former case it is 
dependent on SNR, whereas in the latter case it is not. 

The true envelope detector would yield a Rician distribution (or non-central 


Rayleigh) for the target signal in each of the range cells. Closed from expressions for the 
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probability density function (pdf) of such a distribution exist and involve modified Bessel 
functions [Ref. 8]. The true envelope detector yields a Rayleigh distribution for noise alone. 
However, since the envelope detection approximations x, and X, are being used, the 
derivation of the pdf turns out to be much more difficult. This forms the focus of this 
research. 

The goal is to attempt to determine if a closed form solution can be derived for the 
pdf of the signal produced by (1) and (2). The approach to be taken consists of deriving the 
probability density functions of the smallest components of the equation first. Then the 
transformed density functions have to be derived following the same constructions that 


make up (1) and (2). The pdf of x, is p, (x) , representing the pdf of a range cell with a 
target signal present in it. The pdf of £, is p. (X). It is assumed that the target signal is 


located within the test cell only, and does not contribute power to the noise-only cells 


adjacent to it. 
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HI. PDF FOR THE ENVELOPE DETECTION 
АРРКОХТМАТТОМ аМах{ ШО} + bMin{l,|QI} OF A 
TARGET TEST CELL 


The development outlined in this section derives the density function p.(x) of a 
target test cell sequentially. It begins with the probability density functions of each of the 
I and Q sample expressions. These are known to be Gaussian distributed, being composed 
of known means and normal noise. The densities of the magnitudes of I and Q are found 
next. Then the density functions of the terms Max {|/|,|Q|} and Min {|/|,|Q|} are 
derived. This is followed by the derivation of the density functions of these terms 
multiplied by constants, i.e. aMax {|/|, |Q|} and bMin {|/|,|Q|} . Finally, the probability 
density function of the term x, = aMax {|/|,{Q|} +bMin {\/|,|Q|} is derived. The 
derivation of p, (x) involves a convolution of the two terms in x,. Because the density 
functions of Max and Min operations are quite involved, the convolution of their terms 
generate a very complex probability density function. This may not be the most desirable 
approach. As shown, this density function is such that computational efficiency and 
computer memory capacity become important concerns. Additionally, in order to obtain a 
closed form solution of the convolution integral, curve-fitting techniques are used to model 
the error function. Instead of using what would turn out to be an extremely complex closed 


form solution for p(x), curve-fitting techniques are also investigated to generate 
expressions for p, (x). The expressions for p. (x) are derived for the selections of a and 


b given in Table 1. The coefficients necessary to implement these expressions are tabulated 


to provide a useful design tool for the evaluation of the Pp and to provide a performance 


comparison with other envelope approximations such as £,. 
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A. CLOSED FORM PDF BY ANALYTICAL METHODS 
The inputs into the envelope detection approximation block of Figure 4 are the in- 
phase and quadrature components of the radar return signal. The data sample for each 


channel can be expressed as the sum of a modulated amplitude and some noise, as 
] 2 Acos0 4 nj (9) 


and 


where A is the time-invariant signal amplitude and O is the signal phase. The noise 
variables, nj and no, are assumed to be mutually independent and identically distributed 
Gaussian random variables with zero mean and unity variance (7N(0,1)). By definition, the 
I and Q terms are therefore also normally distributed random variables which are 
independent of each other and identically distributed. As discussed in section II, the signal 
phase is uniformly distributed over the range O...(7x/4). The final result for the 
probability of detection is averaged over the phase. That is, the Pp is calculated as 


ሼ/4| መ= 2 - 
Ро (МК) = : | ЇР, „ (2) ЇР Ne) 9 | | Pc (x, SNR, 8) “< 42 |20 _ (11) 
010 0 Tz 


n 


The average values of / and Q are, respectively, 


т, = 4 050 (12) 
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апа 


то = Аз. (13) 


That is, the noise samples have zero mean. The signal amplitude 1s also related to the signal- 


to-noise ratio (SNR) as A = /202SNR. 
1. PDF of a Gaussian Random Variable 


The probability density functions of / and Q are Gaussian and are expressed as 


[M 
1 26: 
е 


8028.6 4 
с.п (14) 


24 (2) 
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= ፍው] 
1 | 2с? 
е + 


от 


po (q) (15) 


Since the variance is unity, the o2 terms can be omitted in the following derivations 


without any effect. 
2. PDF of the Absolute Value of a Gaussian Random Variable 


Considering only the / random variable, the cumulative distribution function 


(cdf) of / is the integration of the probability density function and can be expressed as 


ዖ,/(0 = Pr(I&i) 
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i -(u-m,) 
P ,(i) = fel zu (16) 


00 


which is the probability that a given realization of the random variable / will not exceed a 
threshold i. The cdf of l/l is then 


Pi, (i) = Pr{{(-) 2 (-i)} 0 [/«ባ) 


= Рг 15!) 


Or 





і (-(и-т)) 
P | (ü) = ie 2 € | 


(17) 


It is known that the probability density function of II is obtained by 
differentiating the cumulative distribution function of Л. Since the integral operation is 
linear, (17) can be expanded into two parts as 

i 


Pn (i) = [pi (u) du 


-і 


ог 
0 і 

P | (i) = [ 2; (и) аи + |p; (u) du . (18) 
0 


- 
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Differentiating this expansion with respect toi, the density function ОҒ ІЛ сап ре expressed 


as 


0 i 
Ри (i) = <р du +£ p, Cu) du (19) 
-i 0 
and becomes 
2 (27) 5 2൮ (2) -7(-2) . (20) 


The final expression for the pdf of IIl is then given by 





-(і-т)) -(-і-т)) 
nos Ll Í Јо 
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መመመ) E DRE 
മഥ ൫ в u (i) (21) 
J2n 
where u(i) is the unit step function. This has to be present to ensure that the variable i, 
representing values of l/l, remains positive in accordance with the absolute value operation. 


In a similar fashion , the probability density function of IQI can be expressed as 


-(q- mg) -(4+ то) 
Ро (4) < | : 24 2 1. (22) 


№ 
= 
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These functions are valid probability density functions since they satisfy the appropriate 
requirements: they are positive for any value of i and q, and they both integrate to unity. 
3. PDF of Max{ll,|Q\} 
To begin, a random variable Z is defined which takes the maximum of two other 


mutually independent random variables X and Y, 


Z = MAX(X,Y) (23) 


The variable Z has a cdf defined in the following manner 


Е,(2) = Рг42 52) 


= Pr (MAX (X, Y) <z) 
Or 


Е7 (2) = Р | Х<г2 7 <г) (24) 


In other words, the cumulative distribution function of Z is the probability that both X < z 


and Y < z for any given value z. Because X and Y are independent from each other, 


Е7(2) = Рг{Х<2} Рг{У< 2} 


or 


Е7 (2) = Ех (2) Ру (2) | (25) 
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` That is, the joint cdf of the two independent random variables X and Y is equal to the 
product of their respective distribution functions, Fy (z) * Fy(z). 

In general, the derivative of a function which is composed of the product of two 
other functions (u and v, for example) which in turn each depend on the same variable (x) 


can be written as [Ref. 9] 


pos zm Du _ (26) 


Applying (26) to (25) to obtain the density function of Z, and knowing that differentiating 
a distribution function with respect to its argument yields the corresponding density 


function, f; (z) 18 obtained as 


S F,(z) = fzx(z) = £ [Fy(z) Fy(2)] 


= Fy (z) 2ғұ(г) +Fy(z) £.Fy(2) 


= Fy (z)fy (2) +Fy(z)fy (z) 
Or 


ይ(።)።= fe(2)Fy(z) +fy(2) Fy(2) | (27) 
This expression represents the pdf of a maximum of two random variables. It can also be 
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derived using [Ref. 8] 


fz (z) = | | f(x, y) dxdy (28) 


where the term AD; is defined as the region of the x-y plane where z < Z < z + dz, and the 


function f(x,y) is the joint density function defined as 


f(x,y) = fy x) fy) (29) 


because of the fact that the two random variables are independent. This approach, however, 
is slightly more complicated. 


In the context of this thesis research, the random variable Z is defined as 


2 = Мах (1, |0]) (30) 


е 


In order to find its pdf using (27), the cumulative distribution functions of the Il and [0] 
random variables are first derived. By definition, these are obtained by integrating their 
respective density functions over the range —© to z. Using (21), the cdf of pjy(i) is given by 
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The function erf is the error function defined as 


х-т т 


^ -г 
Aj | (32) 
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Similarly for IQI, using (22), the cumulative distribution function for Ро is given by 


Po (2) = ര്‌ ae "us ദേ 


Using (27), the density function of Z is then 





pz(z) = Pin (2) Pig (2) +Ро (2) Руд (2) 


or 


ол. 2-2 
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-(2-110) -(2+ то) s 
” : 1, : П 7) erf 2) |u() 


(34) 


That is, the probability density function of the maximum of the two random variables l/l and 
ІОІ 15 a combination of exponentials and error functions. 
4. PDF of Ми, ОП 
The density function of a minimum of two mutually independent random 
variables is found in a similar manner as that for a maximum, discussed in the previous 


subsection. The random variable W is assigned to be 


W = Min(X,Y) (35) 


where X and Y are two independent random variables. The cumulative distribution function 


of W is expressed as 


Fy(w) = Pr{W<w} 


= Pr{MIN(X,Y) sw} 


ОГ 


Fy(w)2 PriX swuYsw] (36) 
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This means that the cdf of the minimum of X and Y is equal to the probability that either 
X < wor Y <w for any given value w of W. Thus the distribution function can be expanded 


as per the definition of an or operator as follows 
Fyw(w) = Pr{Xsw} +Pr{Ysw}-Pr{XswoYsw} . (37) 


Because the random variables X and Y are independent, the cumulative distribution 


function becomes 


Fy (w) = Pr{Xsw} +Pr{Y¥sw} -Pr{Xsw} Pr{Ysw} 


ОГ 


Fy, (w)= Fy(w) +Еу (м) -Ех (№) Ру (№) (38) 


which 15 expressed in terms of the distribution functions of each of the random variables X 
and Y. Differentiating the distribution function of the random variable W with respect to w, 
and thereby applying equation (26) to the third term, the probability density function of W 


can be written as follows 


E Fyw) = fw (w) 


= 4 [Fy (w) + Fy(w) — Ех (№) Ру(») ] 


= < Fx (w) +-«-ዖ› (#) - £ [Fx(w)Fy(w)] 
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= Fy w) & A Fy Qu) - [Fg (ov) Ay (ow) + Руби) Ру (о) | 


= № (№) +Лу (№) - [Fy Ov) fy Ow) + Fy Ow) Fy бе) ] 


= Џубе) -fy (w) Fy (w) 1 + уб») fy w) Fx w)] 
or 
fw (w)= few) П-Ру(») ] +Лу 0) 1 –Рубе)] 39) 
In this thesis it is desired to find the density function of the random variable 
W = Ми (1, 191) (40) 
In accordance with (39), the density function of W is given by 


pw (w) = Pin (w) [1 -Ро (у) ] "Ро (и) ШЕШІ (99) 1 | (41) 


Using (31) and (33) for the cumulative distribution functions of the random variables I/I and 
IQl respectively, and (21) and (22) for the density functions of I/| and ርጋ respectively, the 
probability density function of W can then be expressed as 
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pw (w) = 








-(и-т)) — (и+т,) 
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(ич то)“ - (w+ mo) ü 
ге JM E зе >") 2 jer) 


Or 


pw (w) = 














-(м--т;) -(и+т,) 
ши {Eo : Lacer ജിം 


2/21 12 


5 
I) 





2 


– (» – тр) - (ю+ то) D Я 
ሠ 5 E : 1 (2-60(7- 5) 60 МШ) 2 


(42) 


This defines the probability density function of the minimum of the random variables Л апа 
IQI, expressed in terms of exponentials and error functions. Note that the unit step function 


prevents py,(w) from being negative. 
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5. . PDF of aMax(Ul,Igl and bMin{\N,|Q1} 
Generally, if the random variable Y is a linear transformation of one other 


random variable X of the form 


y = aX (43) 


where a is a constant, the probability density function of Y is then given by [Ref. 8] 


fry) = of) (44) 


а 


where fy is the density function of X. Let the random variable 2" Бе given by 


Z' = aZ = амах (11, Ор. (45) 


The density function p; (z) of Zis given by (34). Assume the constant a is always positive. 


Therefore, using (44), the density function of Z' 1s then 
1 2 
DINE M 


Or 











pz (2) = 2207 
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(46) 


In a similar fashion, let W' be given by 


И” = bW z bMin {1,10} (47) 


where b is a constant. The density function pw (w) of W 15 given by equation (42). Assuming 


that the constant b is always positive, and again using (44), the probability density function 


of W' is given by 


pw (w) = zPw C 


Or 


Ру (w) = 
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(48) 


6. РЕ оғаМах ОВ + РЬМїп{Ш,1О!} 


The envelope detection approximation x, is also arandom variable and is given by 
x, = аМАХ { 11, [9] } «оМШК 111, QI} 


ог 


E (49) 


The density functions of the random variables Z' and W' are known and given by (46) and 
(48), respectively. Since 7 and Q are assumed to be independent random variables, then Z' 
and W' are also independent. This is known because of the fact that at any given time, if Z' 
is a function of l/l, then W' 15 automatically a function of IQI only, and conversely if Z' isa 
function of IQI, then W' is a function of I/I only, by definition of the Max and Min operations 
in xe. 

If two random variables are independent, then the density of their sum equals the 
convolution of their densities [Ref. 8]. Hence the density function of x, is the convolution of 
the density functions Z' and W'. This represents the density function p.(x) of a target test 


cell. The pdf of x, is thus given by 


P(x) = pz (x) @ py, (x) 
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ОГ 


ео 


Ре) = | рр (У) Ри (х= у) ау _ (50) 


This linear convolution integral must be solved to obtain the pdf of xe. Since p> (z) = O0 
for z <0 and pw (w) = O for w <O, then p, (x) = O for x< 0. In addition to removing 


the need for the unit step function in the density function expressions of Z' and W’, the 


integration limits can be set from 0 to x . The density function of x, then becomes 


> (2: > ዘካ) =(2 4m) 
በ ቪር паста 220222 
x ] 2 2 
e 


2-6 х 




















31 


В-н 582 905621 ја | (51) 


As demonstrated, this probability density function 1s a complicated expression. 
An evaluation of this expression is presented next. Each term within the integral is 
expanded. This expansion is necessary to enable a closed form solution to be developed. 


The pdf of x, therefore becomes 
1 2 
p, C) = Sapa | L (C 2EQ – ЗЕН +2ВЈ +285 – ЗЕО +21.0 - 2KA +2BW+ 
0 
KIR+2LM —-2KW+2LN —-2KU -2EG - LDG - BXW - BDW * 2LG — 


LPN + KPS +2LQ -2KV - BDA -2EM + EIO + KIV - LIH + KPU — 


ВР5 +ЕХМ – 110 – 110 + КРЕ +КХЕ – ТРО +КХА +ЕРН + ЕХС + 


2LC —-2KF -2ET—-LDM -LXC * KIU * KXJ - LXT- LPO - LXM - 
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ВІК – ВРК – ІІМ - BIS - EPN - BXA - LPH + KPV + KDA - LDT — 


LXG — BXJ + KIS — BDJ + EXG + EXT -2K/+2LT+ KDW — BXF ~ 


ВРУ ~ ВІО ~ ВРО — BIV + KXW +2BV+2BU+2BR-2EN+2LH + 


KPR +EDT -BDF -2KS -2EC-2KR + EIQ + EPQ + ЕРО + ЕШ + 


ЕрМ-10С-4Кр/4ЕН 4ЕРрС-42В8Е-4ЕРС-42ВА)х 


-1 
зе (-2хиа? + и?Ь? + Зтоа?Б“ +х?а? + и?а? + 2итоа?Ь + 22т2а?Ь + 
а 
е 


2xm,a b * 3mja b  2um;a?b 4 2итаЬ? + 2umoab2 ) ] e 


The variables A...X in equation (52) are symbols which are represented as 


1 
ap (Zumjb * 2um,a 4 moab + umça + xmoa + трађ + ит В) 





А <е 
В f и+ та 
= ег 
| а./2 | 
x (Aum,b + 4хта+ Зтоар + 2итоа + 2хтоа + m;ab + 2итоб) 
С =е 
х-и+тоБ 
Ше ТЕН 
0.2 

= еђј | ——о 

| 2.42 | 
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(52) 


(53) 


(54) 


(55) 


(50) 


(57) 


F = е 


1 2 2 
т (2umob + 2xm,a + moab + umça + xmoa + m) ab + ит) 


2 (4xm,a+ Зтоар +2итоа + 2хтоа + m;ab 2итоБ) 


G 2 е2аБ 


= (2um,b + 2хтоа+ moab + тах+ ита+ m;ab + ит) 


fH =e 
х-и+ тЬ 
а 
1 2 2 
P „аб (2um,a + траб + итоа + тоах + тјаб + итуђ) 





E 


K e erf ( —7 


L = — 


ዐ/2 


пи (4um,b + 4um,a + Зтоа + гитоа + 2хт а + miab+ 2umob) 


M = e2ab 


1 2 2 
Е (2um,b + 2итра + тоаь + тах + ита + таб + ито) 


Nese 
2 = (2итоа+ тра + тах+ита+ та + umob) 
ше 
р f z-u-mjb 
- erf| ————— 
2.2 

1 (2 2 ab rab b 
" aa хтоа + тоаб + тах + ита + тјађ + ито) 


== (4xmça = moab + 2хта + Зтрађ + ита + гитуђ) 
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(58) 


(59) 


(60) 


(61) 


(62) 


(63) 


(64) 


(65) 


(66) 


(67) 


(68) 


(69) 


(70) 


ны (4umoa+ тоар + 2хтуа + Зтјађ + Зита + Зит,ђ) 


шр? (71) 
1 2 2 
LA (um, 3mçab + 2umoa+ 2хтоа + таБ+ 2итоб) 
Tp n (72) 
= (4umob+4xmoa + mgab + 2хта + Зтгаб + 2um,a+2um,b) 
О = е (73) 
= (4итоб + 4umoa + тра +2хта+ 3m;ab + 2um,a+2um,b) 
V = е^ (74) 
А, (2xm,a * траь+ umoa + трах+ m;ab + um,b) 
W = e (75) 
x 16 -и- Mob (76) 
- erf| ————— 
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а. Error Function Approximation 


The density function of the envelope detection approximation x, given by 
(52) cannot be obtained in a mathematical closed form. The specific difficulty in doing this 


is to integrate expressions of the form 


f Cerf (au +b) еси +ди+е) ди (77) 


Or 


? + fu+ 
f {erf (au +b) erf (cu * d) e*" д ŝ}du (78) 
A few closed form solutions exist for indefinite integrals of an error function multiplying 
an exponential, but their variable arguments are simpler and different in form from those 


of (77). Very accurate series expansions offering error residuals to the order of 10° to 
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10-10 also exist for the error function. These expressions, however, are of fifth or higher 
order, and are thus very difficult, if not impossible, to integrate in a closed form. [Ref. 9] 
The simplest approach available to integrate (52) 1n closed form is to 
approximate the error function with a polynomial. A second order polynomial-fitted 
approximation will ensure a closed form solution for the integration of terms in the form 
of (77) and (78). The error function can thus be approximated with a maximum residual 
error magnitude of less than 1.75e-3. In order to achieve this order of residual, the error 
function erf(x) has to be modelled in six segments. For the first five segments the error 


function is approximated by polynomials of the form 


erf (x) чах * bxc (79) 


where the coefficients a, b and c for each segment are outlined in Table 2. For the sixth 


TABLE 2: ERROR FUNCTION APPROXIMATION POLYNOMIAL COEFFICIENTS 


sssi 
1: 0.0-0.5 -0.2539 1.1755 -0.0018 
2: 0.5-1.0 -0.4739 1.3547 -0.0384 


-«3905| 0፡ | 08 | [0 - 


segment, the value of the error function is set to unity, since by constraint the error function 





is very close to a value of one for an argument value beyond 3.0. The comparison of the 
actual error function to its approximations, along with the corresponding residuals, are 
presented in Figures 5 to 14. Note that smaller residuals can be obtained by approximating 
with a higher order of polynomial. À conditional statement is included inside the 


integration loop to verify if the error function argument is negative. This is in order to limit 
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erf(x) 


magnitude of residual 


Ma approximated error function 
actual error function 


0.05 O.1 0.15 0.2 0.25 о.з 0.35 6.4 0.45 


Figure 5: Segment 1 of Error Function Approximation 


0.05 О.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 


Figure 6: Residual of Segment 1 of Error Function Approximation 
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approximated error function 
actual error function 
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Figure 7: Segment 2 of Error Function Approximation 





0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 


Figure 8: Residual of Segment 2 of Error Function Approximation 
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0.96 approximated error function 
actual error function 
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Figure 9: Segment 3 of Error Function Approximation 
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Figure 10: Residual of Segment 3 of Error Function Approximation 
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0.995 approximated error function 
actual error function 
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Figure 11: Segment 4 of Error Function Approximation 
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Figure 12: Residual of Segment 4 of Error Function Approximation 
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Figure 14: Residual of Segment 5 of Error Function Approximation 
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............ approximated error function 


the number of additional integration terms introduced by the error function approximation. 
If the argument is negative, the value of the error function is found from the relation 
erf(—x) = —erf(x) [Ref. 9]. These curves are generated using the program in Appendix 
A. 

As an example of how to integrate (52) using the error function 
approximations, the first term (-2EQ) will be integrated here, for the range where the 
argument of the error function of E (57) is between 0.5 and 1.0 (segment 2). In other words, 


when the integration of (52) (integrating over u) is being performed between the limits 


и = а (то - 1/42) апаи = а (то - J2) , integration of the first term (referred to as 





I1(z)) is expressed as 
а(то- 42) ix (- 22ሀሪ' ተ иЬ + 3moa^b' + 22а? +и?а? + 2итуа 5 + 
I. (z) = | -2ЕОе““ 
1 
а(то- 1/ ‚/2) 


22т ав + 22т ab + 3m;a b + 2um,a°b+2um,ab* + 2um,ab* ) Рт 


2 
-и+ тра -и+ 
а(то- 42) - 4739 ( с J +1354 ( : ጨጨ 


а./2 а,/2 


= --2 е x 


a (mg - 1/ 4/2) 








1 2 2 
ЕЕ (22тоа + тоаб + гта + ита + тј аб + итођ) 


JE (-2гиа? + и?Ь? + Зтоа?Ь? + 2222 + иа? + 2ит„а?Ь + 22тоа?Ь+ 
а 
е 





2:т,аЬ + Зт?а? Б? + 2ита + 2дитав? + 2 b? 
zma m;a ита ита ита ) 1 du 
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or 


0.00005 ; (дгто 2тођ + гт, + тјато – 1 4142та + mrb — 1.4142т,9) 


[1 (2) = | 


(туа + mob) 3 


ра (16848 тор“ 4 16848m; a^ * 5154.435abm, - 9478b^  33696mjamob * 


5754.4355 то ие 


| 0.000025 5; (42то + Атођ + 22т, + 2m,amy — 1.4142m,a+ 2т?Ь - 14142 тов) - 


(туа + тор) 3 


ba (20819т 67 + 20819тга" + 24912.786аБт) - 1895667 + 41638 тјатор * 


24912.7866“то ) | 
(80) 


Using this error function approximation in (52) would result in a closed 


form target test cell density function p, (x) being composed of over 1100 exponential 


terms. Hence, this may be a very costly method of determining p, (x) for a given situation. 


B. CLOSED FORM PDF BY CURVE-FITTING 


The probability density function obtained in Section IILA.6 for the envelope 


detection approximation x, = aMax (I, |Ql) * bMin (|I, |Ql) cannot be practically 


implemented for either system performance analysis or actual operational use. This is due 


to the very complex and tedious procedure required to convert the convolution integral of 


the test cell (52) into a closed form expression. Of course, the integral could be evaluated 


numerically, but this approach would be extremely computationally intensive and time- 
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consuming. This is especially true for the computation of probability of detection curves 

for the GO CFAR processor, which require at least four integrations (see Section II.A). 
The alternative presented here provides a relatively simple closed form expression for 

the density function of the target test cell with arbitrary SNR. Generally, the closed form 


solution for a target test cell density function can be modelled as an exponential of the form 


10 8 7 6 5 4 
p, (x) መ ያሚ *BX «CX e DX e EX Ех + G0 e Hx + 1x +Jx+K (81) 


where the coefficients A, B, ...., K are functions of the SNR and may be expressed as 


polynomials of the form 


A = Ар + As] + Ass നനന 5% А6 (82) 


with similar expansions for B, ...., K. The variable s represents the SNR in decibels. 
Polynomial curve-fitting in a least-squares sense is used to yield (81). Such curve- 
fitting, using a tenth order exponential as shown above, generally provides a maximum 
residual error magnitude of 0.03 if the range of the density function over which the curve- 
fitting is applied is small enough. Now, as is well known, the constraint of “small enough” 


is ambiguous. In this research, the range over which the density function is greater than 


10? was considered the limit of the curve-fitting range. A tenth order exponential is used 
in (81) so as to maintain a minimum residual error. These, however, have the disadvantage 
of not being integrable in closed form, and have to be integrated numerically. Higher order 
exponentials provided no further improvement in the residual error. Smaller order 
exponentials in (81) would ensure that a closed form solution could be obtained for the 


integration of the target test cell density function p.(x) in (7). Such smaller order 


exponentials, however, proved to be very inaccurate. 


The solution of the form (81) is obtained by performing the numerical convolution 
(50) for a number of increasing signal-to-noise ratios and for a given combination of the 
multiplying constants a and b of the envelope detection approximation. Each numerical 
convolution was averaged over the phase angle in the range 0...1/4, as described in 


Section II.B. The results of such an operation are shown in Figure 15, where the 





| 


Figure 15: Test Cell PDF for Increasing SNR 


consecutive density function curves correspond to the test cell density p, for each of the 
signal-to-noise ratios (SNRs) used and where the constants were chosen to be a=1.0 and 
| b=0.25. These particular constant values were chosen arbitrarily in order to demonstrate the 


present discussion. These curves are generated using the program in Appendix B. 
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Each curve of Figure 15 is then curve-fitted, independently of the others, in the form 
of (81). For example, modeling of the first curve (SNR = O dB) over the range 


x = 0.01...5.6 (region where р„ > 1077) yields the curve-fit shown in Figure 16, which 








Actual ре 
Curve-Fitted p, 
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Figure 16: Curve-Fit of p, (x) for SNR = 0 dB, a = 1.0, b = 0.25 
also shows the corresponding magnitude of the residual. The exponential coefficients 


required for this curve-fit are 
A = —0.0015 


B = 0.0459 


-0.5889 


С 


р = 4.262 


Е = –19.1293 
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Е = 55.1303 


С - -102.2169 
H = 119.0197 
- -83.2403 
J = 33.1493 
К = -7.482 


Outside of the curve-fit range, the value of p. (x) is taken to be zero. 

By performing a number of these curve-fits for successive p. (x) curves in Figure 15, 
sequences of coefficients for each of A,B, ....., K can be plotted out. All the curve-fits thus 
obtained are shown in Figure 17, with their corresponding residual error magnitudes. The 
coefficient sequences are plotted out for SNR = 0....21 dB and for a = 1.0 and b 2 0.25 in 
Figures 18 to 28. It 1s rather difficult to extract valid coefficient values from these plots. 
Therefore the exact coefficient values used to produce the plots are tabulated. This is 
actually done for all the envelope detection approximation situations outlined in Table 1 so 
as to provide a good reference design tool. One density p, (x) can be produced at each 
integer value of SNR. Tables 3 to 9 provide all the necessary coefficients. The SNR range 
0....21 dB was chosen since it covers the range over which Py, is less than one (see section 
V). Beyond this range, Pp is constant at one. 

The curves for A,B, ....., K can in turn be curve-fitted 1n the form of (82). For instance, 
the exponential coefficient A of (81) can be expressed as (82) where the coefficients Ag, 
А7,....., Ао аге determined by the polynomial curve-fit. Thus, for a given a and b and SNR 
range, it appears that the density function of a target test cell could be approximated by both 
(81) and (82). This would, in effect, remove the dependence of (81) on the SNR. The 


coefficient curves are also compared to their curve-fits in Figures 18 to 28. As well, these 
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Figure 17: Curve-Fits for p, (x) , with Residual Errors 
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Figure 18: Behaviour of Exponential Coefficient A as a Function of SNR 
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| Figure 19: Behaviour of Exponential Coefficient B as a Function of SNR 
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Residual Error 
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Figure 20: Behaviour of Exponential Coefficient C as a Function of SNR 
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Figure 21: Behaviour of Exponential Coefficient D as a Function of SNR 
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Figure 22: Behaviour of Exponential Coefficient E as a Function of SNR 
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Figure 23: Behaviour of Exponential Coefficient F as a Function of SNR 
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Figure 24: Behaviour of Exponential Coefficient G as a Function of SNR 
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Figure 25: Behaviour of Exponential Coefficient H as a Function of SNR 
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| | Residual Error | 
SNR = 0,1,....., 21 dB Í | 

Actual I | | | 
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Figure 26: Behaviour of Exponential Coefficient I as a Function of SNR 
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Figure 27: Behaviour of Exponential Coefficient J as a Function of SNR 
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Residual Error | 


Figure 28; Behaviour of Exponential Coefficient 85 a Function of SNR 
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show the corresponding residual error magnitudes. Note that the residual error in each case 
hovers very near zero, indicating an effective curve-fit. As it turns out, however, attempts 
to obtain a reasonable target test cell density function using the curve-fits of the coefficients 
A, B.,.....,K were unsuccessful. It appears that the result of a polynomial expansion as in (82) 
is very sensitive to slight changes in its coefficient values. Hence, since a curve-fit 
generally has errors at every point, regardless of how small, the resultant polynomial does 
not behave as expected. 

As mentioned earlier, smaller order exponential expansions approximating p, (x) 
generate significant errors, as can be seen in Figure 29, which shows a third order 
exponential curve-fit. 

The program used to generate all the curve-fits and plots in this discussion is in 


Appendix C. 
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TABLE 3: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 1.0, FOR x, 


Range 
Upper A 
Limit XU 





иии ስ ሽሽ i ስ лии ሽን: т K ലു 


С 

. [2 6.78 -1.9238е-4 | 6.9423е-3 | -1.0802е-1 
| 025 | 70 — -1.2619е-4 | 4.7207е-3 | -7.6157е-2 
ич 7.295 -8.0869е-5 | 3.1468е-3 | -5.2814e-2 
— | 76 -4.76716-5 || 1.93606 -3.3925е-2 


-1.7716е-6 | 8.9455е-5 | -1.9557е-3 
-3.9263е-7 | 2.1412е-5 | -5.0709е-4 


a | om 
0 
0 
7.945 -2.6193е-5 | 1.1146е -2.0473е-2 
1 
2 
7 10.415 -5.3040е-8 | 3.2497е-6 | -8.7064е- 
1 


2 = 
25 - 
3 
-3 
55 | - -3 
2 : - : - 
4 | : : 
05 Е : 
11 
: -6 
-6 


© 


11.105 23166e-8 | -12966e-6 | 3.0154е-5 
11.88 4.6389e-8 | -2.9893e 8.2487e-5 


= 
~ 
~ 
~ 


12.745 4.9775е-8 | -3.5120е-6 1.0648е-4 


Q2 
на 
со 
ኣርን 


4 14.81 4.0119е-8 | -3.4077е-6 | 1.2490е-4 
5 16.035 3.1905е-8 | -2.9665е-6 | 1.1892е-4 


QN 
: 

Сл 

NO 


13.72 4.5677e-8 | -3.5275e-6 1.1728e-4 


5 
5 

17.405 2.1449е-8 | -2.1390е-6 9.0897е-5 
9 


TABLE 3 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b= 1.0, FOR x, 


me је Ге је јо 


1 
2 
3 
5 
7 
10 
11 
12 
13 
14 
15 
| 16 
17 
19 


19 
20 
21 


| 9.3750е+0 | 2.2605е+1 | -3.8683е-9 
aee 


4.7207е-3 
3.1468е-3 
1.9360е-3 
1.1146е-3 
5.8365е-4 
2.6303е-4 
8.9455е-5 
2.1412е-5 
3.2497е-6 
-1.2966е-6 
-2.9893е-6 
-3.5120е-6 
-3.5275е-6 
-3.4077е-6 
-2.9665е-6 
-2.1390е-6 
-1.0544е-6 
3.6074е-9 
8.4043е-7 
1.3038е-6 
1.4038е-6 


-7.6157е-2 
-5.2814е-2 
-3.3925е-2 
-2.0473е-2 
-1.1284е-2 
-5.3839е-3 
-1.9557е-3 
-5.0709е-4 
-8.7064е-5 
3.0154е-5 
8.2487е-5 
1.0648е-4 
1.1728е-4 
1.2490е-4 
1.1892е-4 
9.0897е-5 


4.1129е-5 


-1.8714е-5 





-7.64596-5. 


-1.1726е-4 
-1.3443е-4 


TABLE 3 (cont’d): EXPONENTIAL COEFFICIENTS, а = 1.0, b = 1.0, FOR x, 


| он | азаа | аттен 
эон 
| ањон | зеен | ховын 
പ [ക] വട 
Оз | തര] മ 





































൭൫ 


TABLE 4: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.5, FOR x, 


SNR (dB) 


2 


1 
7 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 

0 

21 


Range 


Upper 
Limit XU 





15485 | 9.112e-10 


15 
15 
02 
03 
55 
155 
3 
37 
55 
44 


Сл 
| i 


Ë 
В 


2.5679е-2 | -3.4668е-1 
1.9372е-2 | -2.7018е-1 
1.2918е-2 | -1.8686е-1 
7.8035е-3 | -1.1764е-1 
4.7479е-3 | -7.4765е-2 
2.7179е-3 | -4.4848е-2 
1.2734е-3 | -2.2174е-2 
4.9103е-4 | -9.0876е-3 
1.2652е-4 | -2.5165е-3 
2.3765е-5 | -5.0945е-4 
5.9043е-6 | -1.2978е-4 
1.6178е-6 | -3.0537е-5 
3.8520е-7 1.3429е-7 
4.5402е-7 | -3.8989е-6 
7.5926е-7 | -1.6489е-5 
6.6902е-7 | -1.5785е-5 
3.0051е-7 | -2.3702е-6 
-3.1705е-7 2.6409е-5 
-1.0975е-6 7.0533е-5 
-1.8827е-6 1.2429е-4 
-2.2461е-6 1.5917е-4 
-1.5160е-6 1.1354е-4 


TABLE 4 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.5, FOR x, 


SNR (dB) 


3.7896е+1 | -7414e«1 | 9.1405е+1 
1 3.2549е+1 | -6.575е+1 | 8.3637е+1 
2.5194е+1 | -5.293е+1 | 7.0114е+1 
1.8047е+1 | -3.970е+1 | 5.5211е+1 


- 
. 
- 


1.3122е+1 | -3.025е+1 | 4.4157е+1 
9.0994е+0 | -2.207е+1 | 3.3975е+1 
5.3428е+0 | -1.381е+1 | 2.2771е+1 
7 2.6724е+0 | -7.451е+0 | 1.3385е+1 
· 9.5011е-1 | -2939е+0 | 5.9730е+0 
2.5533е-1 | -9.0354е-1 | 2.1830е+0 


10 7.3431е-2 | -2.9130е-1 | 8.5036е-1 
11 1.7862е-3 
12 -3.6964е-2 
13 -4.3359е-2 | 2.9147е-1 | -1.150е+0 
14 -3.1541е-2 | 3.0413е-1 | -1.580е+0 
15 -5.1991е-2 | 5.6644е-1 | -3.452640 
16 -1.4429е-1 | 1.4839е+0 | -9.526е+0 


5.5201е-4 
2.0507е-1 


6.7262е-2 
-6.4428е-1 


3 
17 -4.01656-1 || 4.10886+0 | -2.7696+1 
18 -9.7786е-1 | 1.0538е+1 | -7.635е+1 
19 -2.049е+0 | 2.3915е+1 | -1.893е+2 
20 -3.339е+0 | 4.2635е+1 | -3.7096+2 
21 -2.818е+0 


3.8142е+1 | -3.524е+2 


TABLE 4 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.5, FOR x, 


roman [+ |) | 
о [эт жжет | ээп 
2.7182е+1 | -8.458е+0 
2.4447е+1 | -8.708е+0 
. 4. | -056/541| 2.2285е+1 | -9.119е+0 
1.9990е+1 | -9.631е+0 
| 6. | 2360e4 | 1.6866е+1 | -1.016е+1 
1.36746+1 || -1.0836+1 
| 8 _ | -8343e40 | 1.0434е+1 | -1.164е+1 
_ 9 | 409760 | 8.3157е+0 | -1.289е+1 
7.4870е+0 | -1.4766+1 
6.80706+0 | -1.7096+1 
5.5877е+0 | -1964е+1 





3.7456е+0 | -2.230е+1 
-2.4275е-2 | -2386e«1 
-1.600е+1 | -1.509е+1 
-7,4596+1 
-2.8116+2 
-9.5736+2 
-2.9196+3 
-7.009e43 


3.7069е+1 
2.6044е+2 
1.0984е+3 
3.8492е+3 
1.0380е+4 


-7.587е+3 1.2057е+4 
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TABLE 5: EXPONENTIAL COEFFICIENTS, a= 1.0, > < 0.25, FOR x, 


Range 


Upper 
Limit XU 


a 
С) 





__ о | o | 56 | -1.5422е-3 | 4.5899е-2 | -5.8885e-1 
еми _ -1.1524е-3 | 3.5342е-2 | -4.6724е-1 
. 05 | 596 | -7.3691е-4 | 2.3410е-2 | -3.2073e-1 

| 40 | 61 — -4.6829е-4 | 1.5451е-2 | -2.1993е-1 
| 025 | 6.5 -2.9259е-4 | 1.0057е-2 | -1.4916е-1 
| (0 | 615 - -1.5489е-4 | 5.5750е-3 | -8.6645е-2 
| 065 | 7095 -7.4596е-5 | 2.8221е-3 | -4.6143е-2 
፳ከ፡ በኡ ፳ -3.0750е-5 | 1.2297е-3 | -2.1284е-2 
| uo ሎክ | -8.6030е-6 | 3.6777е-4 | -6.8269е-3 

| 45 | 8:21 | -1.3960е-6 | 6.4832е-5 | -1.3188е-3 

а MES 6.4795е-8 | -2.7352е-6 | 4.0778е-5 
2.3171е-7 | -1.2180е-5 | 2.7321е-4 
2.4897е-7 | -1.4384е-5 | 3.5862е-4 
2.2852е-7 | -1.4434е-5 | 3.9556е-4 
1.7959е-7 | -1.2307е-5 | 3.6627е-4 
1.0541е-7 | -7.5467е-6 | 2.3071е-4 
7.7617е-9 | 5.1223е-7 | -7.3435е-5 
-1.0992е-7 | 1.2479е-5 | -6.2665е-4 
-2.3077е-7 | 2.7493е-5 | -1.4599е-3 

-3.3934е-7 | 4.4222е-5 | -2.5740е-3 

-4.1982е-7 | 6.0504е-5 | -3.8995е-3 
-4.6637е-7 | 7.4716е-5 | -5.3582е-3 


TABLE 5 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.25, FOR x, 


we] 9 | « [ ፡ | 5 [.።. 
Г 0 | 42620е+0 
3.485 1e+0 4.7874е+1 1.09676+2 
2.4803е+0 3.6714е+1 9.11116+1 
1.76786+0 2.83176+1 7.6320е+1 


1.2496е+0 | -6.5036+0 2.1771е+1 | -4.701е-1 6.3951е+1 


7.61356-1 || -4.1606+0 || 1.46496+1 4.79306+1 
4.2708е-1 9.16526+0 3.3833е+1 
2.0893е-1 5.0835е+0 2.1641е+1 
202075 50 2.0726е+0 1.0776е+1 
1.5444е-2 5.6590е-1 3.9729е+0 
-1.4745е-4 4.2063е-2 8.70026-1 
-3.3979е-3 -1.1715е-1 -4.71016-1 
-5.0431e-3 -2.4549e-1 -2.052e40 

l 


0 
1 
2 


aà | р-а 


м 
(ታን 


-6.1599е-3 -3.8036e-1 -4.316е+0 
-6.1978е-3 -4.51016-1 -5.873е+0 
-3.8945е-3 -2.2216е-1 1.0567е+0 
3.0544е-3 9.21406-1 4.5690е+1 
1.8324е-2 4.3796е+0 2.1960е+2 
4.5471е-2 | -9.1937е-1 | 1.2599е+1 7.5399е--2 
8.808 1е-2 2.9692e+1 2.1882е+3 
1.4796е-1 6.1622е+1 5.6515е+3 
2.2647е-1 1.1743е+2 1.3479е+4 


> 


— 


5 
6 
7 


кез 


— 





-1.962е+0 


E 


9 
0 


പ്ര 


ኩጋ 
— 


jà 
оо 


ርጋ 
UJ 


TABLE 5 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.25, FOR x, 





- 
















эле 
элем 


TABLE 6: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 3/8, FOR x, 





- 
С) 


-1.0551е-3 | 3.2194е-2 | -4.2355е-1 
-7.8236е-4 | 2.4621е-2 | -3.3409е-1 
-5.0488е-4 | 1.6466е-2 | -2.3164е-1 
-3.2244е-4 | 1.0926е-2 | -1.5975е-1 
-1.8673е-4 | 6.6024е-3 | -1.0077е-1 


3.7489е-3 | -5.9940е-2 
1.8189е-3 | -3.0660е-2 
7.4499е-4 | -1.3316e-2 
1.9928е-4 | -3.8292е-3 
3.4402е-5 | -7.1993е-4 
2.9103е-6 | -6.3560е-5 
-4.2046е-6 1.0972е-4 
-6.5844е-6 1.8570е-4 
-7.6542е-6 | 2.3624е-4 
-8.2503е-6 | 2.8031е-4 
-8.6974е-6 | 3.2652е-4 
-8.7357е-6 3.6289е-4 
-7.9157е-6 3.6312е-4 
-5.8704е-6 2.9418е-4 


1 


015 
015 
02 
025 
035 
-4.6655е-5 
135 -1.8049е-5 
265 -4.5077е-6 


7 


10 -6.0718е-8 
12 1.0133е-7 
13 1.0839е-7 
14 1.0667e-7 
15 1.0218е-7 
16 9.3045е-8 
17 | 6310 | 1485 | 7.6490е-8 
EET: 5.1897е-8 
19 2.1102е-8 | -2.4552е-6 | 1.2384е-4 | 
0 10.505 -1.2266е-8 | 2.18816-6 | -1.7090е-4 | 
21 12:255: -4.2599e-8 | 7.4224е-6 | -5.7898е-4 


. 


TABLE 6 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 3/8, FOR x, 


2.5587е+0 | -1.216е+1 3.7145е+1 | -7.312е+1 9.0702е+1 


X 
| 1 - 


1 
2 
3 
12 -2.9748е-3 2.9754е-2 | -1.9258е- 8.1016е-1 -2.164е+0 
2 
1 
1 


а тож 
5 E 


TABLE 6 (cont’d): EXPONENTIAL COEFFICIENTS, a= 1.0, b = 3/8, FOR x, 


ЕСЕ БЕГЕН БЕГЕН БЕГЕН 
Го | Patient | зове | Tato 
നി്‌ 


+ 
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TABLE 7: EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8, FOR x, 


Range Range 
SNR (dB) 


Lower Upper 
Limit Х|, Limit XU 





- 
e 


-1.5822e-3 -5,9773е-1 
-1.0277e-3 -4.1472e-1 
-6.6134e-4 -2.8683e-1 
-4.2092e-4 -1.9723e-1 
-2.6300e-4 -1.3393e-1 
-1.3964е-4 -7.8090е-2 
-6.3592е-5 -3.9441е-2 
. 13 
1665 1 
| 2525 _ 


0 
0 0971е- 
0 = 
| “-2.39686-5 -1.6722е-2 
-6.0010е-6 -4.8188е-3 
| -9.5220e-7 -8.9789е-4 
8 -9.327e-8 -9.1379e-5 
1.33 7.3921е-8 1.1780e-4 
1.885 1.1898e-7 2.1161е-4 
2.525 1.2892е-7 2.7120e-4 
1.2802e-7 3.2398e-4 
1.2580e-7 3.8669e-4 
1.1921e-7 4.4761e-4 
1.0417e-7 4.7828e-4 
7.9951e-8 | -89591e-6. 4.4593е-4 
4.7268е-8 3.0943е-4 
9.6952e-9 3.6398e-5 


11.905 20.46 -2.6865е-8 4.7589е-6 | -3.7602е-4 


7 
10 
11 
12 
13 
14 
15 
17 
18 
19 
20 
21 
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TABLE 7 (cont’d): EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8, FOR x, 


3.2776е+1 | -6.524е+1 | 8.1971е+1 
2.5382е+1 | -5.264е+1 | 6.8951е+1 
1.9563е+1 | -4.235е+1 | 5.7924е+1 
1.3220е+1 | -3.016е+1 | 4.3596е+1 
7.8944е+0 | -1.916е+1 | 2.9609е+1 
7 4.0704е+0 | -1.064е+1 | 1.7857е+1 
1.5092е+0 | -4.373е+0 | 8.2716е+0 
3.9300е-1 | -1.325е+0 | 2.9972е+0 
10 6.5609е-2 | -2.9716е-1 | 9.2877е-1 
11 -7.9292е-2 | 2.5466е-1 | -4.2467е-1 
12 -2.0899е-1 | 8.6341е-1 | -2.259е+0 
13 -3.7300е-1 | 1.7824е+0 | -5.5576+0 
14 -6.2056е-1 | 3.3838е+0 | -1.2216+1 
15 -1.035е+0 | 6.4082е+0 | -2.6516+1 
16 -1.669е+0 | 1.1671е+1 | -5.484е+1 
17 -2.462е+0 | 1.9297е+1 | -1.020е+2 
18 -3.098е+0 | 2.6888е+1 | -1.574е+2 
19 1.9560е-1 | -2.643е+0 | 2.4160е+1 | -1.468е+2 
20 8.84426-1 | -1.648е+1 | 1.8126е+2 
21 1.0822е+1 | -1.530е+2 | 1.4694е+3 


TABLE 7 (cont’d): EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8, FOR x, 


неј ој s 
БЕНШЕРГЛЕРІЛЕГІІ 
ടു പര 
ക്കി മ asas 
а | азан овен | этен 
ване 
6 ажиыа | елен വമ 
Оз [эзет [эзет силен 
പ [വ വോ 
9 ടട മോ 
10 [ase Tome | ied 
La ве നാ 
se |e = 
6:28 
m m 
е 
— 
ш 
NN 
N 
Шог 








ШЕЛІ 
1221 
IIT 










NA 


-5.530е+4 
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TABLE 8: EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393, FOR x, 












Range 
Lower 
Limit XL 


Range 
Upper 
Limit XU 






SNR (dB) 


Б 
В 


-1.8377е-3 | 5.3571е-2 | -6.7308е-1 
5.665 | -1.1871е-3 | 3.5775е-2 | -4.6488е-1 
-7.5975е-4 | 2.3744е-2 | -3.2007е-1 
-4.8463е-4 | 1.5736е-2 | -2.2044е-1 
-3.0099е-4 | 1.0193е-2 | -1.4896е-1 
-1.5977е-4 | 5.6676е-3 | -8.6809е-2 
-7.2123е-5 | 2.6945е-3 | -4.3514е-2 
-2.7062е-5 | 1.0714е-3 | -1.8365е-2 
-6.7088е-6 | 2.8441е-4 | -5.2386е-3 
-1.0933е-6 | 4.9927е-5 | -9.9742е-4 
-1.4511е-7 | 6.6944е-6 | -1.3671е-4 
| ass | 9148 -— 4.3950e-8 | -3.0324e-6 | 8.5833е-5 
12 9.5767e-8 | -6.3203e-6 | 1.8027е-4 
13 1.0725е-7 | -7.5820е-6 | 2.3393е-4 
1.0874е-7 | -8.3681е-6 | 2.8277е-4 
15 1.1220е-7 | -9.4587е-6 | 3.5192е-4 
1.1575е-7 | -1.0742е-5 | 4.4164е-4 


17 14.255 1.1395е-7 | -1.1682е-5 | 5.3215е-4 
15.5 1.0529е-7 


-1.1953е-5 6.0411е-4 
8.9424е-8 | -1.1243е-5 | 6.3014е-4 
20 6.6838е-8 | -9.2765е-6 | 5.7415е-4 
21 3.9626e-8 


-5.9821e-6 4.0160e-4 
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TABLE 8 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393, FOR x, 


3.9790e- -2.267е+0 8.3527е+0 | -1.999е+1 3.0472е+1 
7 1.7864е- -1.086е+0 4.2861е+0 | -1.105е+1 1.8296е+1 
5.5014е-2 | -3.6341е-1 | .1.5717е+0 | -4.492е-0 8.3925e+0 


1 
1 
1 
узе 


21 -1.5768е-2 4.0030е-1 -6.853е+0 7.9908е+1 | -6.246е+2 


ሆሪ. 


TABLE 8 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393, FOR x, 


we o 
о [тыш [эзш Эзе 
зво ыа | ляна 
2. | -640261 | 287666+1 | -8205е+0 
. 3. |2864 | 26850641 | 48550560 
.- 4. | -:964| 251146+1| -9.040е0 
. 5. | 45% | 222076+1| -95086+0 
Cs 1 ጨ[መ men 
E -1.060е+1 -1.152е+1 
[ക IET] 41272641 
= 
OB | 
പ 
അ 
ഈത്‌ 
ഇത്‌ 
അ 
ഈത്‌ 
അ 
ഇന്ന്‌ 








ламен 
вен 
ТІЛІ 


1.0771е+4 
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TABLE 9: EXPONENTIAL COEFFICIENTS, a = 0.96043, b = 0.39782, FOR x, 


у= ы 
-1.4228е-3 -5.3667е- 
-1.0458е-3 -4.2043е-1 
| -6.1273е-4 | 2.1292е-2 | -2.9066е-1 
6.165 | -4.2793е-4 -1.9977е-1 
-2.4809e-4 -1.2606e-1 


133464 | 4796403 | 7445352 


Range Range 
Lower Upper 
Limit xy, Limit xij 


SNR (dB) 


C 
1 





-6.4488е-5 | 2.4384е-3 | -3.9852е-2 
2 3286е-5 | 9.3387е-4 | -1.6217e2 


-5.8691е-6 | 2.5196е-4 | -4.6999е-3 
-9.4667е-7 | 4.3803е-5 | -8.8678е-4 


o 


3 
3 
4 
5 
870 -1.3013е-7 | 6.0822е-6 | -1.2580е-4 
11 1.335 9.295 3.9229е-8 | -2.7342е-6 | 7.8260е-5 
12 9.925 8.4119е-8 | -5.6251е-6 | 16255е-4 
6 
6 
6 
6 
5 
-5 


13 9.4973е-8 | -6.7975е-6 | 2.1234е-4 
14 9.6305е-8 2.5675е-4 
15 9.9680e-8 3.2058e-4 
16 1.0205e-7 3.9933e-4 
17 1.0037e-7 4.8078e-4 
9.2926e-8 5.4675e-4 
19 7.8618e-8 5.6808e-4 
10.18 5.8443e-8 5.1454e-4 


l IL87 |] 20475 4 3417262 17552174e:6 | 315415000 


| 
| 
: 


| 
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TABLE 9 (cont’d): EXPONENTIAL COEFFICIENTS, а = 0.96043, b = 0.39782, FOR x, 


о | вов тузеп [эзеп [элее [тосо] 


3.6878е-1 -2.127е+0 7.9267е+0 -1.919е+1 2.9595е-1 
-9.8490е-1 | 3.9402е+0 1.7310e+1 


4.9987е-2 | -3.3447е-1 1.4654е+0 | -4.244e40 8.0344е+0 
1.0341е-2 | -7.6957е-2 3.8215е-1 -1.284е+0 2.9082е+0 


1 ട്‌ 
-2.6576e-3 | 2.7033&2 | -1.7723е- -2.002e40 
à 2 

















TI .- 

EE | 

E — 

9.5835е-! 

3.40276-1 

| | 

4.919040 

-1.0890+1 

-2.4656+1 

5.4896+1 | 

2.14656+1 
17 -1.6396е-2 | 3.1887e-1 | -4.200е+0 -2.3216+2 

је-2 1" 409232: -6 6.0547е+1 | -4.140е42, 






-1.8919е-2 


-1.8922е-2 4.5206е-1 -7.327е+0 8.1574е+1 | -6.157е+2 
-1.4054е-2 3.6043е-1 -6.229е-0 7.3261е+1 | -5.768е+2 


— 
NO 


20 
21 
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TABLE 9 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.96043, b = 0.39782 FOR x, 





тэ -7.579е+1 | 3.1419е+1 | -7.680е+0 
-7.158e+1 | 3.0683е+1 | -7.988е+0 










Се [шашы [ен [чег 
С он | tees | sen 
СУ зөввө | sse | immer 
[9 







1.8372e43 | -4.770е+3 | 5.46186+3 
3.0134е+3 | -8.625е+3 | 1.0908е+4 
2.8923е+3 | -8.259е+3 | 9.9716ሬ+3 
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Figure 29: Third Order Curve-Fits for p(x), with Residual Errors 
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IV. PDF FOR THE ENVELOPE DETECTION 
APPROXIMATION alll + bIQ| OF A TARGET TEST 
CELL 


A. CLOSED FORM PDF BY ANALYTICAL METHODS 

The envelope detection approximation £, given by (2) is a simplified version of x,, 
the approximation considered so far in this thesis. It does not have the Max and Min 
operations and can therefore be calculated faster than x,. Analytically, X, also has the 
advantage that the probability density function can be derived without the need for error 
functions. As can be recalled in the pdf development for x,, the error function terms arose 
from the expressions for the cumulative distribution functions. These were required to 
describe the statistical distribution of the Max and Min operations. 

In this situation, the derivation of the density function of X, for a target test cell 
requires the convolution of two simple linear transformations of the random variables I/I 
and IQI by the multiplying coefficients a and b. Applying (44) to equations (21) and (22), 
the two terms in the convolution are 


1 ድ 2 1 ድ 2 
(8) 1 pum Шато @ 
р = 





апа 


2 2, Бе 
Рио പ്പം 25 С (84) 


These can be convolved since / and Q are assumed to be mutually independent Gaussian 
random variables. As can be recalled from Section II, the convolution arises from the 


addition of the two independent components (ІЛ апа БО! 
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The convolution integral is given by 


p,(%) = | pan (y) py o 2-У) ду 


--60 


or 


~ 


x 


p.) = [Рал (у)Рьо 2-У) ду. (85) 
0 


The integration limits reduce as shown when the unit step functions of (83) and (84) are 


removed from the integration integral. Expanding (85), the convolution integral becomes 


Ен 10+ 2ут 4 E: ፡| F ~ \ 
е a фес ? -+е- б 


ЈХ 





А ( P... + Е 
е 


2 1 
Р, Е 22ዕ፳| | ም у. 
0 


2 
- (- 2Хуа” = yb + m;a b? + a+ ya +2 ymga’b + moa b! + 2хт oa b + 2ут jab’) 


1 
то 
е 245 dy 





ОГ 





x : [- (а? + b^) y+ (4m,ab’ + 242 - 28860 ዕ - 2m,ab’) y+ (47 т а? Ь - 
р, (8) = | је“ 
с 2abn | 


2 
m;a b -Ха - ജുമ -2Ятра Ь) | Tae 





Í 2 22 [- (2 + 2) у + (Ата? + 4тра?Ь + 229 - 2тра?в- 2трав“) у+ (- та?» - 
а 

fe 

0 


$ d - mid b –2атрађ ) ] 20 


fi») 


| [— (а? + b) y + (2ха? - 2mga b - 2m,ab’) y+ (A&mga b-m;a b. - 


2:5 
e^? 





СУ У യയ >> 


ያ d - mia b 2 атрађ ) ] dy А 





1 
22 
era b 


0 


[- (а +?) у? + (4тоаБ' + 23а – 2тра’Ь -2т а) у+ (- m;a b! - 


„2 А 
Х а”-түа б-23туа b ) lay | | (86) 


Each of the four integral terms in (86) are now of the form 





ደ 
[മ (87) 
0 
which integrate in closed form to 
4CA - 
de l——— B 242 +ይ 
[ a сво er (88) 


2/-4. 2./-А 2 БА 


Applied to (86), this gives the result 


28 1 Jn аа В, 2AR +B, 4 
Ре С Завп| 2 „ГА “| 23) АЕ =й 














‚т ха ел В, ЛЭЭ! И 
2 J-A 2 J-A 


4ር,4 - 8 И 
Jn Мы. j мыл) ё 


2/-4. 216) “2 
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where 





22 - 2 
T 4A 


А = 


T Ва / 242 + В, 
ДІ seen == 
2:74 | 24—À 








DO 2 2 1 
В, - (4moab +2£a – 2тоа ђ— 2таб) (Е рога) 


= (4&mga^b - т/а?» -ха 


25 252 
C, = (= туар 


С; = (4Ята? в-та? р? -Ха -moa b -2imga^ Б) 3 


Са = (- mab 


„22 


#а?- -moa 2, затоа“ 0) ( 


) 


2476 


22222 Oe О 


2. а? -moa 2,2 -25тоа” Ж ORO? 
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B, (4m,ab! 4 242 - 2тоа?Ь = 2таЬ*) ж 
2a b 


1 
22055 





a^ p? 


2) 


2,2 


В, - (Атјађ“ + 4тоа 5 + 2ла? - 2тоа b - 2m,ab’) 5) 


) 


QUI Rude. 2 1 
-moa b - 247702 ത്തു; 
2а Ь 


(89) 


(90) 


(91) 


(92) 


(93) 


(94) 


(95) 


(96) 


(97) 


(98) 


В. CLOSED FORM PDF BY CURVE-FITTING 
Using the same procedure as that used for Section III.B, the density function p, (X) 


is also obtained by numerical convolution for anumber of signal-to-noise ratios. The results 
are curve-fitted using an exponential expansion of the form of (81). The exact coefficient 


values for each SNR are in Tables 10 to 16 in order to provide a quick and easy reference. 
This procedure allows the two envelope approximations x, and x, to be compared in terms 


of their representing coefficients. 
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TABLE 10: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 1.0 FOR £, 


Range 
Upper A 
Limit XU 


-5.6415е-5 | 2.2342е-3 | -3.8108е-2 
-4.2441е-5 | 1.7303е-3 | -3.0380е-2 
-3.1021е-5 | 1.3057е-3 | -2 3666е-2 

8 


റ 


SNR (dB) 





-1.3325e-2 


HELM 

| 

О» [вз јаз взен 
0 
‚2 
8 
2 


51 6.9408е-4 
| 025 9.785 -2.5321е-6 -2.9351е-3 
09. 10.305 | -6.6428е-7 | 3.6566е-5 | -8.6903е-4 
BEEN 10895 | .23983e-8| 14881e-6 | -4.1031e-5 


3 

4 

4 
1.425 12.305 2.740Зе-8 -1.8370е-6 5.2550е-5 
13.145 2.7140е-8 -1.9782е-6 6.1884е-5 

7 

6 


5 
05 
3 
05 
5 
05 
1 
25 
45 
35 
09 


02 28 11.560 2 2540е-8 | -13770е-6 | 3.5360е-5 


| 285 — 
-7.95956-9 | 1.6944е-6 | -1.6077е-4 


8.915 -9.9361e-6 4.6842e- -9.5072e-3 


2.835 14.090 2.3047e-8 | -1.8092e-6 6.0928e-5 


2 


1 
7 
10 
11 
[2 
13 
14 
15 
16 
17 
18 
19 

0 
7! 
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TABLE 10 (cont'd): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 1.0 FOR &, 


SNR (dB) ан ЕЕ. н 

-1.841е+0 2.8377e+1 
-1.528е+0 2.5808е+1 
аа ea 
-9.9219e-1 2.0573e+1 
-7.6809e-1 1.7904е+1 
-5.1141е-1 1.3703е+1 
ТІСІ ков 
-9.7447e-2 3.8972e40 
4.0194e-3 -6.1962е-2 
8.1102е-3 -4.9075е-1 
1.1741е-2 -1.072е+0 
1.3298е-2 -1.318е+0 
-5.7953е-3 1.1571e41 
-9.0311e-2 1.2689e42 
-1.5848e-1 2.9271е+2 
19 -2.3215е-1 | 4.6394е+0 | -6.327e+1 | 5.8138е+2. 
-3.2342е-1 1.5195е+3 


84” 


TABLE 10 (cont’d): EXPONENTIAL COEFFICIENTS, a= 1.0, b = 1.0 FOR x, 


ет [о к 

541419 
5.713640 
:6093e40 
“| даден | 12915041 | -6576е40 
-7.1896+0 
. 61 | -152:1| [1:10926+1| -785060 
-8.633e+0 
8 | 53960 | 6816250 | 9493640 
. 9. | -1290 | 472406+0| -106%; 
-1.259041 
4.1940е40 | -1.494е+1 
-ህ — 


зен 
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TABLE 11: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.5 FOR £, 


Range 
Upper A C 
Limit xy 


005 -8.6196е-2 





5 
Š 


.005 


з 
M] 
UJ 


-1.1347е-7 6.4902е-6 | -1.6294е-4 


= = 


e 
л 
: 
A 
= 


2 -5.5374e-7 | 2.7465e-5 | -5.9376e-4 
1.14 10.01 -9.2300e 5.4751е-6 | -1.4238е-4 


12 1.775 10.68 -1.1528е-8 7.3671е-7 | -2.0688е-5 
13 2.505 11.435 6.037е-1 -1.5928е-8 8.5076е-7 
е 


-8 

1 

- 7 
-9 6 
-9 
10 


20 -8.8275е-7 | 5.2265е-5 
21 3.4981е-7 | -4.3519е-5 


TABLE 11 (cont'd): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.5 FOR x, 


mel s | s | r | 1 2 
Го Геи ome) rier [amet | Seer 
-2.60e+1 | 3.70206+1 
ET 

2.7514е-3 4.7734е+2 


1 
7 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 


TABLE II (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.5 FOR x, 


2 
4 
5 








элее 
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TABLE 12: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.25 FOR x, 


Range 


SNR (dB) Upper 


Limit Хү) 





- 


005 -8.88 16e-4 
005 -6.6626е-4 
005 -4.8777е-4 
005 -3.4597е-4 
005 -2.36016-4 
005 -1.55956-4 
005 -9.8414е-5 
01 -5.0651е-5 
025 -2.1057е-5 
105 -5.4285е-6 
375 -3.3928е-6 
795 -4.7802е-7 
1.38 4.4594e-8 
2.04 -1.521e-10 
2.715 -1.0876e-9 
3.6 -2.5077e-9 
4.525 -5.6096е-9 
5.55 -1.2042е-8 
67 | -2.441е-8 
7.985 -4.6017e-8 
9.425 -7.9046e-8 
11.03 -1.2136e-7 


© 


IA 


99 


~ | СА 


ኩጋ 


0 


~ 


1.6444е-4 


-3.4782е-1 
-2.7696е-1 
-2.1623е-1 
-1.6450е-1 
-1.2116е-1 
-8.6829е-2 
-5.985 le-2 
-3.4111е-2 
-1.5975е-2 
-4.9015е-3 
-3.4603е-3 
-5.6183е-4 

6.8575е-5 
4.0578е-8 
-1.9784е-6 
-6.3354е-6 
-1.8818е-5 
-5.2465е-5 
-1.3639е-4 
-3.2740е-4 
-7.1324е-4 


1.9393е-5 | -1.3860е-3 


TABLE 12 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, р = 0.25 РОК x, 


ај е Ге Гр" Ге Ги |. 
_- |» | [=] 
| 

1 


1.6358е-1 -1.042е+0 | 4.2776е40 | -1.134е-41 1.8905е+1 
5.5182е-2 | -3.8912е-1 1.7836е+0 | -5.336е+0 


1.0174е+1 


3.2817е- -1.959e+0 7.5072е+0 | -1.849е+1 2.8422е+1 


Е ЁОЛ 
Lm 
| 4 
-2.1818е-2 8.6865е+0 
2| m 


; - 1 
3 - 
-4.7141е-6 1.0104е- -1.0702е-3 6.0659е-3 | -1.8244е-2 
5 
4 - 
3 : - 
3.9680е-3 | -7.4579е- 9.45 19е-1 -8.175е+0 4.7653е+1 
2 ን 


1.08086-2 | -2.31516-1 3.3602е+0 | -3.345е+1 2.2555е+2 
2.66 13e- -6.4615e-1 1.0663е+1 | -1.211е+2 9.3400е+ 
5.8329е-2 -1.600е+0 2.99006+1 | -3.8526+2 3.3789е+3 


кә 


1 
7 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 


TABLE 12 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 0.25 FOR х, 


е јо ро 
Ls Ген тен ames 
ste 
ШЕНИЕ ЛЕССЛЕТ 
а эзоњи | inen | «amen 
т [шен [сыш osea 
— 
Се 
ЯЛ. 
се 
EE 
Се в 
Lu |ጋ 
Сва 
Се Де 
ШЕШЕ: 
Са В 





-3.4700e+0 -1.305е+1 
ах 
эзин 
122 


-9.759е+4 
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TABLE 13: EXPONENTIAL COEFFICIENTS, a = 1.0, b = 3/8 FOR x, 


Range 
Upper A 
Limit XU 





-2.4855e- 


5.865 -5.9627е-4 | 1.8544e-2 
| 60 — -4.4671е-4 | 1.4312е-2 | -19764е- 
71 


: 


С 
1 


-3.2481е- 1.0754е- -1.5346е-1 
-2.3131е 7.9286е-3 | -1.1713е-1 


-1.5888е-4 | 5.6558е-3 | -8.6757е-2 


| «5 | 69 | -1.0527е-4 | 3.9019e-3 | -62319е-2 


2 


7.665 | 73.41466-51 1389563 | 271011 


8.515 -112366-6| 546916-5 | TRS 


© 


-6.6350е- 2.5708е- -4.2912е- 


Фр 
= ==> 27 
Ou. 
wo 
ሠ 
à 
= 
о 
== 


8.065 -1.1382е-5 | 4.9132е-4 | -9.1576е-3 


3 


10.23 4.6154е-9 | -3.0650е-7 | 8.9859е-6 


3 
3 
3 

-3.0786е-7 | 1.7138e-5 | -4.1691е-4 
7 
9 


14 3.105 11.76 -2.835е- 1.9028е- -6.2129е-8 


17 14.84 3.854е- -5.8870e-9 | 3.4002е- 


4 
-4 
4 
5 
5 
505 -6.5642e-7 | 3.4302е-5 | -7.8293е-4 
13 2.315 10.95 1.1272е-9 | -7.7951е-8 | 2.3756e-6 
1 
1 
1 
1 | 
0 
0 


8356-1 

15 12.67 -4.250е-11 | 2.9915е-9 | -9.9091е-8 

16 4.985 13.69 | .3.001൦-11 | 1.7139е-9 | -3.9367е-8 
1 


7 
6 


17.575 6.477е- -8.8982e-8 | 5.4011е-6 
7 


21 12.06 21.03 3.0399е-9 | -5.0390е-7 3.7284е-5 


18 7.345 16.125 2.262e-1 -2.938 le-8 1.654 1e- 
1 


TABLE 13 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, b = 3/8 FOR &, 


(бше 


1.0246е-2 | -8.4819е-2 4.6215е-1 -1.670е+0 3.9248е+0 


! 


5.8185е-3 | -5.1414е-2 2.9937е-1 -1.158е+0 2.9180е+0 
-1.5279е-4 1.6637е-3 | -1.2098е-2 5.9686е-2 | -1.9986е-1 


1 
эе 


-4.1764е-5 4.6513е-4 | -3.4065е-3 1.6807е-2 | -5.8636е-2 


1 
7 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 


93 


TABLE 13 (cont’d): EXPONENTIAL COEFFICIENTS, a = 1.0, 6 = 3/8 FOR x, 


-3.749е+1 1.6903е+1 -5.320е+0 





ሚሙ” 
4 | 2፡6:| (ጠመህ| 588 
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TABLE 14: EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8 FOR X, 


Range 
Upper 


Limit XU 


790: 
E жие 
അ me 
— 
| 





2.3358е-2 -3.0497е-1 


-4.4683e-5 -3.0179е-2 
-1.4787e-5 -1.1256e-2 
-1.4578e-6 -1.4159e-3 
8.78 -7.4567e-7 -8.5150е-4 
-3.89776-7 -5.0221e-4 
12 | 15 | 996 1.359e-10 1.7391e-6 
3 | 2255 | 1066 | 1.5245е-9 3.0647е-6 


-1.539е-10 -3.4041е-7 
15 -1.856е-10 -4.8170е-7 
-1.652е-10 -4.7352е-7 
17 -3.901е-11 6.6338е-8 
3.066e-10 2.1770e-6 

1.0559е-9 8.2474е-6 
20 2.5021е-9 2.3456е-5 
21 5.0595е-9 | -8.1189е-7 | 5.8164е-5 


TABLE 14 (cont’d): EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8 FOR x, 


е је Ге ја 


2.2460e+0 -5.678е+1 | 6.7297е+1 
1.8517е+0 | -8.703е+0 | 2.6268е+1 | -5.100е+1 | 6.2110е+1 


1.485 1е+0 | -7.213е-0 2.2486е+1 | -4.507е+1 5.6570е+1 
1.1725е+0 | -5.894е-0 1.9009е+1 | -3.939е+1 5.1028е+1 
9.0511е-1 -4.720е+0 1.5783е+1 | -3.388е+1 4.5409е+1 


3.96376+1 
3.3911е+1 
2.4961е+1 
1.35 10е+1 
4.2050е+0 


6.7543е-1 | -3.667е+0 | 1.2759е+1 | -2.848е+1 
4.87056-1 | -2761е+0 | 1.003281 | -2.336е+1 
2.9035е-1 | -1.731e+0 | 6.62186+0 | -1.626ሌ+1 
-7.3540е-1 |`° 3.02256+0 | -8.0376+0 
1.6699е-2 | -1.2450e-1 | 6.1069e-1 | -1.986е+0 
1.0914е-2 | -8.8563e-2 | 4.7347е-1 | -1.680е+0 | 3.8842е40 
6.8405е-3 | -5.9020е 3.3577е-1 | -1.270е+0 | 3.1320е+0 
-5 
4 


9020е-2 

2.6969е-2 | -1.0335е-1 
2.0040е-2 | -6.9377е-2 
-6.8804е-3 | 2.6616е-2 
-1.2917е-2 | 6.0961е-2 
-1.9178е-2 | 1.0926е-1 
| -34005e-6 | 5.7572е-5 | -2.7059e-4 -3.7666е-3 | 5.8466е-2 
1.4813е-1 | -8.3369е-1 
9.1942е-1 | -6.226е+0 
4.1260е+0 | -3.191е+1 
1.6043е+1 | -1.400e42 


1 
7 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
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TABLE 14 (cont’d): EXPONENTIAL COEFFICIENTS, a = 31/32, b = 3/8 FOR £, 





n 


3.8734e43 





oF 


TABLE 15: EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393 FOR 7, 


Range 


SNR (dB) Lower 


Limit XI, 


> 
e) 





5.645 -8744е-4 | 2.6152е-2| -3.370е-1 
-6.586е-4 | 2.0281е -2.6911e-1 


-3.405e-4 -1.593e-1 
-2.343e-4 -1.182е-1 
| 5 | 675 -1.553e-4 -8.493e-2 
| «5 | 7025 | -9.921e-5 -5.913e-2 
BER. И -1.691е-5 -1.252е- 
|^ 455 | 819 -1.514e-6 -1421e-3 


-6.080e-7 3.1075e -6.9511e-4 


8.675 


s Ja ЧЕНЕ 
ፍስ | ሂስ Сл | ©л | ол | л | Ул | tA л 


‹ 
a | = 


54 т -4.807е-4 | 1.5289е-2 | -2.096е-1 


-1168e-8 | 6.6877е -1.669e-5 
13 2.245 10.53 1.5193е-9 | -1.036 3 
11.305 -4.70е-10 | 3.1455е-8 |  -9321e-7 


15 -4.17е-10 -9.330е-7 
16 -1.71е-10 -1.88е-7 
17 3.798е-10 2.3450е-6 
18 1.4311е-9 8.9045е-6 


16.895 3.2174е-9 | -4.135е-7 | 236138 


2 
3 
E 
2 
Е 
-3.832е-7 | 2.0648е-5 | -4867е-4 
w ር 
7 
8 
-9 


5 


20 18.46 5.9566е-9 | -8.427е-7 | 53146e-5 


21 11.615 20.215 9.7702е-9 -1.530е-6 1.0702е- 
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TABLE 15 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393 FOR £, 


н 6 | ዩ | ፣ [ 9 [ 8 - 


Го | ама бет | ут өтт 
7 
4 
4 


1 

-2.7081е- 5.2928е-3 -6.9543е-2 6.2408е-1 
| 1 
{ 1 


4 
-7.8867е- 1.7058е-2 -2.4976е- 2.5099e+0 
-1.9670е-3 4.7312е-2 -7.7288е- 8.6887е+0 
-4.4011е-3 1.178 4е-1 -2.147е+0 2.6953е+1 


1 
2 
3 
5 
7 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
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6.9432е+1 


6.4009 
5.8369е+1 
5.259бе+1 
4.67636+1 
4.08726+1 
3.5082е+1 
2.5757е+1 
1.3855е+1 
3.9184е+0 
3.2817е+0 
2.9915e+0 
1.1409e-1 
-7.6762e-2 
3.1192e-2 
5.0992е-2 
-1.6839е-2 
-6.3258е-1 
-3.797е+0 
-1.712е+1 
-6.643e+1 
-2.304e+2 


+ 


Ф 
кеі 


TABLE 15 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.948, b = 0.393 FOR X, 


на со : 


-4.8616+1 1.90376+1 -4.423e40 









- 
© 





-4 595e+1 1.8558e+1 -4.660е 


-4 001е+1 1.7459e+1 -5.342е+0 


-3.676е+1 1.6850е+1 -5.827е+0 


-4.308е+1 1.8030е+1 | -4.960е+0 
-3.334е+1 1.6207е+1 | -6.444е+0 





-2.9836+1 -7.227е+0 
-2.3376+1 -8.087е+0 
-1418e«1 -9.024e40 
-5.591e40 -1.018е+1 
-5.255640 -1.235е+1 
-5.080е+0 -1.497е+1 
-6.2182е-1 -1.664е+1 
-2.8564е-1 -2.0386+1 
-5.1125е-1 -2.537е+1 
-5.5837e-1 -3.1436+1 
-3.47636-1 -3.858е+1 
1.8330е+0 | 4.3987е+0 | -4.363е+1 
1.46096+1 -2326e41 
7.5616е+1 1.5776e42 
3.3041e42 1.1792е+3 
1.2834е+3 6.0691е+3 


0 
1 
2 
3 


м | = | pat 





— 


5 
6 
17 
8 
9 
20 
21 


= 


мә | m 
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TABLE 16: EXPONENTIAL COEFFICIENTS, a = 0.96043 р = 0.39783 FOR x, 






Range 
Upper 
Limit XU 





-7.7426e-4 


51 -5.4233е-7 


-3.3557е-7 
-9.5417е-9 


13 2.28 1.3047e-9 
14 -4.066е-10 
15 -3.630е-10 
16 -1.531е-10 
17 3.247е-10 
18 1.2409e-9 
19 2.7948е-9 
20 5.2091е-9 
21 8.5565е-9 


2.3443е-2 
1.8099е-2 
1.3692е-2 
1.0085е-2 
7.1837е-3 
4.9808е-3 
3.2954е-3 
1.7854е-3 
6.3191е-4 
6.4184е-5 
2.8055е-5 
1.83 19е-5 
5.5151е-7 
-9.0119е-8 
2.7606е-8 
2.6193е-8 
8.9322е-9 
-3.9778е-8 


-6.3508е 


ВИН! 
| б 
ഡി: 
CA ы 
ж 
о) 
4> 


2.0683е 


-3.6393е-7 
-7.4668е-7 
-1.3578е-6 


9.6223е 
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-3.0587е-1 
-2.4325е-1 
-1.9001e-1 
-1.4490e-1 
-1.0721e-1 
-7.7361e-2 
-5.3479e-2 
-3.0428e-2 
-1.1431e-2 
-1.3129e-3 


-4 


-4.3746e-4 
-1.3887e-5 

2.7462e-6 
-8.2906e-7 
-8.3554е-7 
-1.8131е-7 


-6 


-1.4935е-7 7.9300е-6 
2.1058е-5 
4.7713е-5 


-5 


TABLE 16 (cont’d): EXPONENTIAL COEFFICIENTS, a = 0.96043 b = 0.39783 FOR 
Хе 


кај 010551 0562 
Го | о тако эки Заден | яаг 
1.00466+1 3.3813е+1 

-1.14006-1 3.8132е+0 
-1.344е+0 | 3.18056+0 


13 -4.8227е-5 5.3596е-4 | -3.9345е-3 1.9834е-2 | -7.2238е-2 


2 
Зак 
21 -4.0094е-3 -2213e42 
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TABLE 16 (cont’d): EXPONENTIAL EE a = 0.96043 b = 0.39783 FOR 


=e е 
Го | виа | вне 
Е “4.674040 
МЕ -4.208е+1 -4.975e40 
9 
шэг 





-4,4386+0 





















[eme 
4 1.6660е+1 | -5.843е+0 
s [imas | ier | rane 
ടി sume 
СИ ЕЕЕ 
അന്നാ 
О ame юше обе 
[ക ine 
ie [ase a sin 
ወ [መባ ኤ። = 
നന നി 


6.0649e+3 
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V. COMPARISON OF THE APPROXIMATIONS 
амах ПЛО) + МіО) AND alll + bIQI OF A 
TARGET TEST CELL 


Signals detected by the envelope detection approximations x, and £, given by (1) and 
(2) have significantly different probability density functions. This fact is true for both cases 


of noise and noise plus target. À graphical comparison of x, and X, with the density 


function of the true envelope detector x = „12 + Q? is shown in Figure 30 for the case 





Figure 30: PDF for x, x, and X,, No Target, n = 1 


in which no target is present. These curves are generated by numerical convolution using 


the program in Appendix D. The density of x, is a result of the convolution of (46) and (48). 
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The density of x, is obtained by convolving equations (83) and (84). Finally, the density 
function of x is obtained from [Ref. 8] 


2 


Eco 
Р ሠ” - хе “и(х) . (99) 


The coded sections for each of these are identified in the program. The multiplying 
coefficients of both approximations are arbitrarily set at a — 1.0 and b — 0.25. 

The results of the numerical convolutions shown in Figure 30 are verified against 
experimental results. Data sets generated for / and Q are sorted into histograms, shown in 


Figure 31. The star symbols “ * " represent the true envelope detector x, while the thinner 


О 
(0 







° 
N 


[у 


1 


Figure 31: Histograms of x, x, and x,, No Target, n = 1 


© 
ന 


ር 
о 


Number of Values in Each Segment, Normalized 
© © 
A 


ከ | Ш 
Value of Trials 


bars represent x, and the wider bars represent X,. The program in Appendix E is used to 


generate these histograms. 
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As shown in Figure 30, the envelope detector approximation x, models Ше шие 


envelope detector the best. The approximation £, is significantly shifted from the true 
detector. These curves are used to compute the function D in (7) and (8), for the Pp and 
Ше Ред respectively. As explained in Section ILC, p, , is an n-fold convolution of these 
curves, for a given n reference cells. 


The disparity between £, and x for the function Ру, п is enhanced as n is increased. 


This can be seen in Figures 32 ad 33, where each of x, x, and x, has been convolved 


0.3 


0.25 


0.1 


0.05 





Figure 32: PDF for x, x, and <,, No Target, n = 4 


for n = 4 and n = 16. These n-fold convolutions are generated using the program at 


Appendix F. The approximation x, generally has the same mean as the true detector x. The 


mean of the approximation £,, however, is increasingly offset from that of the true detector 
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as n increases. For n = 1, the mean of X, is offset from x by about 0.5 dB in SNR. For n = 


4, this offset is about 1.25 dB, whereas for n = 16, the offset is approximately 4 dB. 


0.08r - 





x 


Figure 33: PDF for x, x, and £,, No Target, n - 16 


Such differences ultimately create significant differences in the Pp and the Ред СО 
CFAR performance, as compared to the values obtained for a true envelope detector. As an 
example, the experimental Pp curves for each of x, x, and £, are plotted in Figure 34. These 


are generated using the Monte Carlo simulation program at Appendix G for a = 1.0 and b 
= 0.25 and n = 4. This results of this simulation program are validated against similar 


simulations performed in previous research [Ref. 3, Ref. 4]. The thresholds for each 
detector are chosen so as to yield a Pra = 10^. The threshold for Х, 18 selected from 


available plots [Ref. 1] as T = 5.537. The thresholds for x and x, are determined from 
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experimental trials using the program at Appendix H. They are set at 7, = 3.5 and 


T, = 4.2. The Pp curves show that Хх, provides the most degraded detection sensitivity 


0.9 
0.8 


0.7 


probability of detection 
о о о 
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Figure 34: Probability of Detection Curves for x, x, and X,, a = 1.0, b = 0.25, n 4. 


compared to the true envelope detector. The approximation x, provides a significantly 
better Pp performance in comparison. As n is increased, this advantage becomes even more 


pronounced. 


108 


VI. CONCLUSION 


A. SUMMARY 

The main contribution of this thesis is to determine if the probability of detection for 
а СО CFAR radar processor using the envelope detection approximation 
х, = аМах (|1,|0|) + bMin {|/|,|Q|} can be expressed analytically in a closed form. 
The solution to this problem depends solely on whether the probability density function of 
a radar range cell p, can be derived in a closed form for both of the cases where a target 15 
present and where no target is present. It is shown in Section III that this cannot be done 
without the use of either an error function approximation or an exponential curve-fitting of 
the density function. 

The closed form solution based on an error function approximation is very complex. 
As stated in беспоп Ш.А.6.а, it is comprised of more than 1100 first order exponential 
terms. For the closed form evaluation of the Pp, each of these terms needs to be integrated 
in closed form. This approach is not practical. Therefore, the best option is to evaluate the 
Pp numerically once the error function approximation has been applied to p... Although this 
is still very complex and tedious to program, it has the advantage of being one of the most 


accurate methods of computing the Pp analytically, as long as the residual error in the error 
function approximation is kept to a minimum. Another advantage is that this approach is 
valid for any value of SNR and multiplying coefficients a and b. It can therefore also be 
used to compute the Pra. 

The exponential curve-fitting approach involves the curve-fitting of a number of p, 
curves generated numerically for a given range of signal-to-noise ratios. The exponential 
coefficients resulting from these curve-fits are then summarized in reference tables. This 


has the advantage of enabling a rapid implementation of p, within the Pp expression. Then 
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the Pp can be evaluated numerically, as per the error function approximation approach. The 
main disadvantage of this method is its lack of versatility. It is dependent upon the 
multiplying coefficients a and b of the envelope approximation. Furthermore it can only be 
solved for whole number values of SNR (i.e. 0, 1, 2,.....), up to the SNR required for the 
Рр = 1. 

An additional envelope detection approximation £, = a|/| + b|Q| is considered in 


Section IV. It is offered as an alternative to x,. Since it is a simplified version of x,, its 


density function p, (X) is expressed relatively simply in a closed form solution. It is also 
expressed in the form of exponential curve-fits (as done for x,), with the required 
coefficients summarized in reference tables. 

Finally, the two envelope detection approximations x, and x, are compared to each 


other and to the true envelope detector x. This is done in terms of the density functions for 


non-target range cells and in terms of the Pp performance. In both instances, it is shown 
that the approximation x, delivers a performance which is significantly better than that 
provided by $,. Ultimately, the choice of envelope detection approximation to be used is 
to be based on Pp performance and the ease of implementation. Although x, performs 


better in terms of detection capabilities, X, is simpler to implement since the Max and Min 


operations do not have to be performed. 


B. FUTURE RESEARCH RECOMMENDATIONS 
The possibilities for future research on GO CFAR performance using the envelope 
detection approximation x, = aMax {{/|,|Q|} +bMin {|/|,|Q|} are quite attractive from 


a mathematical point of view. The ultimate goal of such research would be to provide the 


radar design and analysis community with a simple method to compute the Pp and the Ppa 


theoretically. This solution would have to be independent of the signal-to-noise ratio in the 


test cell (or reference cell) and independent of the multiplying coefficients a and b. 
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Simplified analytical solutions for the Pp and the Pra are invariably contingent upon 
the ability to integrate expressions of the form of equations (77) and (78) in closed form. 
This presents a challenging mathematical problem. Solving it would certainly prove to be 
very useful to anyone using this particular approximation x, for an envelope detector, 
whether it is fora GO CFAR processor or for some other use. 

Although less glamorous in its possibilities, the ability to approximate the error 
function in a single segment —© to © while maintaining a very small residual error would 
be very useful. This would greatly reduce the complexity of the closed form probability 
density function of a test cell as presented in Section IIL6. This could be done by 
appropriately modulating either an inverse tangent function or a hyperbolic tangent 
function. À closed form solution would then be required for the integral of this function 
multiplied by an exponential of second order. 

Finally, the generation of the Pp curves using the analytical solutions presented in this 
thesis would be very useful. This could be done for both types of envelope detection 
approximations x, and £,, for all the cases listed in Table 1. The number of reference cells 
used should be n = 1, 2, 4, 8, 16 and 32. These represent the majority of the situations 


encountered in GO CFAR processing. 
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APPENDIX A - MATLAB PROGRAM FOR 
POLYNOMIAL CURVE-FITTING OF ERROR 
FUNCTION 


%polynomial curve-fitting of error function 


xz.01:.01:.5; 
y-erf(x); 
pl-polyfit(x,y,2); 
f=polyval(p1,x); 
e=abs(f-y); 
plot(x,y,x,f,’:”), grid 
xlabel(‘x’) 
ylabel(‘erf(x)’) 


Zosubplot(212) 
Zoplot(x,e),grid 
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APPENDIX B - C PROGRAM FOR NUMERICAL 
CONVOLUTION FOR p, OVER A RANGE OF SNR 


/*probability density functions of various envelope detectors 
foaGOCFAR */ 


#include <stdio.h> 
#include <math.h> 
#define pi (4.0*atan(1.0)) 


main(){ 
FILE *Pcell, *PcellArea; 
int idelay,c,cl,c2,x,dataLength,s; 
float n; 
double pdfCellTgt[16000], pdfSample,pMax[ 16000], 
pMin[16000],phi,sum, 
u,snr,ss,a,b,m1,m2,e1,e2,e3,e4,ef1,ef2,ef3,ef4; 


dataLength- 16000; 

pdfSample-1.0/200.0; 

aero: 

6=0.25;. 

if ((Pcell=fopen(“Pyn2.dat”,”w”))==NULL){ 
printf(“????”); 
exit(1); ] 

for (s=0;s<=29;s++){ 


snr=pow(10.0,((float)s/10.0)); 
ss-sqrt(2*snr); 


for (c=0;c<dataLength;c++) pdfCellTgt[c}=0.0; 
for (phi=0;phi<=9*p1i/40;phi+=pi/40){ 
m1=ss*cos(phi); 


m2=ss*sin(phi); 


for (c=0;c<dataLength;c++){ 
u=(double)c*pdfSample; 


[*envelope detection approximation aMax (IIl,IQI ) -bMin (IIIl,IQI) */ 
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el=exp(-pow((u/a-m 1),2.0)/2.0)+exp(-pow((u/a+m 1),2.0)/2.0); 
e2=exp(-pow((u/a-m2),2.0)/2.0)+exp(-pow((u/a+m2),2.0)/2.0); 
e3=exp(-pow((u/b-m 1),2.0)/2.0)+exp(-pow((u/b+m 1),2.0)/2.0); 
e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
ef 1=erf((u/a-m2)/sqrt(2.0))+erf((u/a+m2)/sqrt(2.0)); 
ef2=erf((u/a-m1)/sqrt(2.0))+erf((u/a+m 1)/sqrt(2.0)); 
ef3=2.0-erf((u/b-m2)/sqrt(2.0))-erf((u/b+m2)/sqrt(2.0)); 
ef4=2.0-erf((u/b-m1)/sqrt(2.0))-erf((u/b+m 1)/sqrt(2.0)); 
pMax{[c]=(e1*ef1 + e2*ef2)/(2.0*a*sqrt(2.0*pi)); 
pMin[c]=(e3*ef3 + e4*ef4)/(2.0*b*sqrt(2.0*pi)); 


/*envelope detection approximation alll + b[QI */ 


ы el=exp(-pow((u/a-m 1),2.0)/2.0)+exp(-pow((u/a+m1),2.0)/2.0); 
e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
pMin[c]=e4/(b*sqrt(2.0*pi)); 
pMax[c]=e1/(a*sqrt(2.0*pi)); */ 


/*envelope detector sqrt(IIIA2 + IQI42) */ 


js pdfCellTgt[c]-u*exp(-pow(u,2.0)/2); */ 
} 


idelay=0; 
for (cz0;c«dataLength;4---c) ( 
idelay+=1; 
sum=0.0; 
for (cl=1;cl<=idelay;++c1) 
sum+=(pMax([cl]*pMin[idelay-c1+1]); 
pdfCellTgt[c]+=(sum*pdfSample); 
} 


} 18 ОЕ PHASE ООРТ Ж+Ж ЕН 


for (x=0;x<dataLength;x++){ 
pdfCellT gt[xJ=pdfCellT gt[x]/10.0; 
fprintf(Pcell,”%e\n” pdfCellTgt[x]); 
} 


} /***** END OF SNR Т.ООРт”тжжжжжэжжэ жээ жи ээн эээ 
fclose(Pcell); 


} /ፆዎችትችችችች ЖЕ ЕРЕ ЖЕЕ РЕ ЕЧ ЕЕ жжжжжжжж END MAIN ж/ 
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APPENDIX C - MATLAB PROGRAM FOR 
EXPONENTIAL CURVE-FITTING OF p, FOR A 
RANGE OF SNR 


Yocurve-fitting of target test cell pdf curves over 
доа range of SNR values (in dB) 


load el 
с=е 1; 
clg 
clear p 
hold 
order=10; 
SNR=21; 
sampleSize=1/200; 
for q=1:SNR+1, 
г1(9)=1; 
while c(r1(q),q)<.00001, 
r1(q)=r1(q)+1; 
end 
г2(4)<г1(4)+100; 
while c(r2(q),q)>.00001, 
r2(q)=r2(q)+1; 
end 
x=[r1(q)-1:r2(q)-1]’.*sampleSize; 
p(q,:)=polyfit(x,log(c(r1(q):r2(q),q)),order); 
f=exp(polyval(p(q,:),x)); | 
plot(x,c(r 1 (9):г2(9) а), E X,abs(f-c(r1 (9):г2(9),9))) 
епа 
hold 
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APPENDIX D - C PROGRAM FOR p, WITH NO 


TARGET 


/*probability density functions of various envelope detectors 
for a GO CFAR for no target case */ 


. #include <stdio.h> 


#include <math.h> 
#define pi (4.0*atan(1.0)) 


main(){ 
FILE *Pcell, *PcellArea; 
int idelay,c,c1,c2,x,dataLength,s; 
float n; 
double pdfCellTgt[16000] ችን. ыш. 
pMin[16000],phi,sum, 
u,snr,ss,a,b,m1,m2,e1,e2,e3,e4,ef1,ef2, ef3, ef4; 


dataLength- 16000; 

pdfSample=1.0/200.0; 

a=1.0; 

b=0.25; 

if ((Pcell=fopen(“Pyn2.dat’,”’w”))==NULL){ 
рги (“7777”); 
ехш1); ) 

ss=0.0; 

for (c=0;c<dataLength;c++) pdfCellTgt[c]=0.0; 


ml=ss; 
m2=ss; 


for (c=0;c<dataLength;c++){ 
u=(double)c*pdfSample; 


/*envelope detection approximation aMax (III,IOI --bMin (III,IQI) */ 


e 1=exp(-pow((u/a-m 1),2.0)/2.0)+exp(-pow((u/a+m 1),2.0)/2.0); 
e2=exp(-pow((u/a-m2),2.0)/2.0)+exp(-pow((u/a+m2),2.0)/2.0); 
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e3=exp(-pow((u/b-m 1),2.0)/2.0)+exp(-pow((u/b+m 1),2.0)/2.0); 
e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
ef l=erf((u/a-m2)/sqrt(2.0))+erf((u/a+m2)/sqrt(2.0)); 
ef2-erf((u/a-m1)/sqrt(2.0))--erf((u/a4-m 1)/sqrt(2.0)); 
ef3=2.0-erf((u/b-m2)/sqrt(2.0))-erf((u/b+m2)/sqrt(2.0)); 
ef4=2.0-erf((u/b-m1)/sqrt(2.0))-erf((u/b+m 1)/sqrt(2.0)); 
pMax[c]=(el *ef1 + e2*ef2)/(2.0*a*sqrt(2.0*p1)); 
pMin[c]=(e3*ef3 + e4*ef4)/(2.0*b*sqrt(2.0*p1)); 


/*envelope detection approximation alll + bIQ! */ 


/* | elzexp(-pow((u/a-m1),2.0)/2.0)*exp(-pow((u/a4-m1),2.0)/2.0); 
e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
pMin[c]=e4/(b*sqrt(2.0*p1)); 
pMax[c]=e1/(a*sqrt(2.0*p1)); */ 


/*envelope detector sqrt(IIl^2 ч-10172) */ 


/*  pdfCellTgt[c]-u*expCpow(u,2.0/2); */ 
} 


idelay=0; 
for (c=0;c<dataLength;++c) { 
idelay+=1; 
sum=0.0; 
for (cl=1;cl<=idelay;++cl) 
sum+=(pMax[cl]*pMin[idelay-cl+1)]); 
pdfCellTgt[c]+=(sum*pdfSample); 
} 


for (x=0;x<dataLength;x++) { 
pdfCellT gt[xJ=pdfCellTgt[x]/10.0; 
fprintf(Pcell,”%e\n” pdfCellTgt[x]); 
} 


fclose(Pcell); 


} ከህ. ЕК ЖА ЕЕ ЕЖ RE EER EEE EERE EERE AEE EERE EERE k END MAIN ВИ 


117 


APPENDIX E - MATLAB PROGRAM COMPARISON 
OF ENVELOPE DETECTION APPROXIMATION 
METHODS BY HISTOGRAM 


Yocomparison of various envelope detectors 
gZofor GO CFAR by histogram 


clear 

clg 

6=.25; 

а=1; 

М=30000; 
q-randn(1,N); 
i=randn(1,N); 
Q-abs(q); 

I-abs(1); 

z-—sqrt(1.^2 -- q.^2); 
zi-a*max(I,Q)2b*min(I,Q); 
z2-a*] -- b*Q; 
[m,n]-hist(z,200); 
[m1,n1]-hist(z1,200); 
[m2,n2]-hist(z2,100); 
bar(n,m./((1/48)*N),’*’) 
hold 

bar(n1,m1 /((1/48)*N)) 
bar(n2,m2./((1/24)*N)) 
hold 
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APPENDIX F - C PROGRAM FOR z-FOLD 
CONVOLUTION FOR p, , 


/*n-fold convolution of reference cell PDF for a GO-CFAR */ 


#include <stdio.h> 
#include <math.h> 
#define pi (4.0*atan(1.0)) 


main(){ 


ү, 


FILE *Pcell, *PcellArea; 

int idelay,c,c1,c2,x,dataLength; 

float s,n; 

double pdfCellTgt[16000],pdfSample,phaseSample,pMax[ 16000], 
p2[16000],p1[16000],p3[16000],p4[ 16000],p5 [106000], 
pMin[16000],phi,theta,sumPhase[16000],area,sum, 
u,snr,ss,a,b,m1,m2,e1,e2,e3,e4,ef1,ef2,ef3,ef4, 


dataLength- 16000; 
pdfSample=1.0/200.0; 
a=1.0; 

b=.25; 


if ((Pcell=fopen(“pyn.dat”,”’w”))==NULL){ 
printf(“??7?”); 
ех (1); ) 


ss=0.0; 
በ11=55:; 
m2=ss; 


for (c=0;c<dataLength;c++){ 
u=(double)c*pdfSample; 


/*envelope detection approximation aMax(IIl,IOl) -bMin (IIl,IQI) */ 


el=exp(-pow((u/a-m1),2.0)/2.0)+exp(-pow((u/a+m1),2.0)/2.0); 
e2=exp(-pow((u/a-m2),2.0)/2.0)+exp(-pow((u/a+m2),2.0)/2.0); 
e3=exp(-pow((u/b-m 1),2.0)/2.0)+exp(-pow((u/b+m 1),2.0)/2.0); 
e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
ef 1=erf((u/a-m2)/sqrt(2.0))+erf((u/a+m2)/sqrt(2.0)); 
ef2=erf((u/a-m1)/syrt(2.0))+erf((u/a+m 1)/sqrt(2.0)); 
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ef3=2.0-erf((u/b-m2)/sqrt(2.0))-erf((u/b+m2)/sqrt(2.0)); 
ef4=2.0-erf((u/b-m 1)/sqrt(2.0))-erf((u/b+m 1)/sqrt(2.0)); 
pMax(c]=(el *ef1 + e2*ef2)/(2.0*a*sqrt(2.0*pi)); 

pMin[(c]=(e3*ef3 + e4*ef4)/(2.0*b*sqrt(2.0*p1)); — */ 


/*envelope detection approximation alll + bIQI */ 


/*  el=exp(-pow((u/a-m1),2.0)/2.0)+exp(-pow((u/a+m1),2.0)/2.0); 
e4=exp(-pow((u/b-m2),2.0)/2.0)+exp(-pow((u/b+m2),2.0)/2.0); 
pMin[c}]=e4/(b*sqrt(2.0*pi)); 
pMax[c]=e1/(a*sqrt(2.0*p1)); */ 


/*true envelope detector sqrt(IIIA2 + IQIA2) */ 
pl[c]=u*exp(-pow(u,2.0)/2); 


} 


/*do not perform this convolution for true envelope detector */ 
|“ idelay=0; 
for (c20;c«dataLength;--4c)( 
рїЇсЇ 0.0, 
idelay+=1; 
sum=0.0; 
for (cl=1;cl<=idelay;++c1) 
sum+=(pMax([cl]*pMin[idelay-c1+1]); 
pi[c]=sum*pdfSample; 
} 7. 


/*convolve up to here for n=1 */ 


idelay=0; 

for (c=0;c<dataLength;++c){ 
р2[сј=0.0; 
idelay+=1; 
sum=0.0; 
for (cl=1;cl<=idelay;++c1) 

ѕит+=(р1[с1]*р1(еіау-с1+1]); 

p2[c]-sum*pdfSample; 

} 


/*convolve up to here for n=2 */ 
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idelay=0; 
for (c=0;c<dataLength;++c) { 


p3[c]=0.0; 

idelay+=1; 

sum=0.0; 

for (cl=1;cl<=idelay;++c1) 
sum+=(p2[cl]*p2[idelay-cl+1]); 

p3[c]=sum*pdfSample; 


/*convolve up to here for n=4 */ 


/* 


idelay=0; 
for (c=0;c<dataLength;++c){ 


p4[c]=0.0; 
idelay+=1; 
sum=0.0; 
for (cl=1;cl<=idelay;++c1) 
sum+=(p3[cl]*p3[idelay-c1+1)); 
p4[c]=sum*pdfSample; 
{рип  (Рсе|,”%е\п”,р4[с]); */ 


/*convolve up to here forn=8 */ 


idelay=0; 
for (c=0;c<dataLength;++<c) { 


p5[c]=0.0; 

idelay+=1; 

sum=0.0; 

for (cl=1;cl<=idelay;++c1) 
sum+=(p4[cl]*p4[idelay-c1+1]); 

p5[c]=sum*pdfSample; 

fprintf(Pcell,”%e\n”,p5[c]); 


/*convolve up to here for nz16 */ 


fclose(Pcell); 


] മാട END MAIN */ 
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APPENDIX G - MATLAB PROGRAM FOR MONTE 
CARLO SIMULATIONS OF Py 


%GO CFAR Monte Carlo simulation for probability of detection 
clear 


Т=3.5: 

11:42: 

T2=5.537; 

ብ! 

b=.25; 

n=4; 

S=10000; 
PHI-pi/4; 
phi=0:PHI/20:PHI; 
зпг=0:25; 


w-2*n; 


for B=1:length(phi), 
for C=1:length(snr), 
SNR=104(snr(C)/10); 
A-sqrt(2* SNR); 
d=0; 
41=0; 
d2=0; 
X-zeros(1,5); 


I-A*cos(phi(B)) + randn(1,S); 
Q=A*sin(phi(B)) + randn(1,S); 
=sqrt(I.A2 + Q.42); 
X l=a*max(abs(I),abs(Q)) + b*min(abs(I),abs(Q)); 
X2=a* abs(I) + b*abs(Q); 


i-randn(S,w), 

q-randn(S,w); 

x—sqrt(1.^2 + 1.72); 

x l=a*max(abs(1),abs(q)) + b*min(abs(i),abs(q)); 
x2=a*abs(i) + b*abs(q); 


Id=sum(x(:,1:n)’); 


lg=sum(x(:,n+1:w)’); 
z-max(ld,lg); 
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Vt-T*z/n; 


ld12sum(x1(:,1:n)"); 
lg1=sum(x1(:,n+1:w)’); 
z1=max(ld1,lg1); 
Vtl=T1*zI1/n; 


142:55011(х2(:,1:0) 
]g2=sum(x2(:,n+1:w)'); 
z2-max(1d2,1g2); 
Vt2=T2*z2/n; 


d=length(find((X-Vt)>0)); 
Pd(B,C)=d/S; 
dl-zlength(find((X1-Vt1)20)); 
Pd1(B,C)=d1/S; 
d2=length(find((X2-Vt2)>0)); 
Pd2(B,C)=d2/S; 
end 
end 


PD=mean(Pd); 
PD 1=mean(Fd1); 
PD2=mean(Pd2); 


save PD PD PD1 PD2 
quit 
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APPENDIX H - MATLAB PROGRAM FOR MONTE 
CARLO SIMULATIONS OF Pra 


GO CFAR Monte Carlo simulation for probability of 
false alarm 
clear 


qug 
11542: 
алі; 
0=529. 
п=4; 
S=100000; 


\=2*п: 


d=0; 
91=0; 
X=zeros(1,S); 


І-тапап(1,5); 

Q=randn(1,S); 

X=sqrt(L.A2 4 Q.^2); 

X1=a*max(abs(1),abs(Q)) + b*min(abs(1),abs(Q)); 


i-randn(S,w); 

q-randn(S,w); 

x-sqrt(1.^2 + 1.^2); 

x 1=a*max(abs(i),abs(q)) + b*min(abs(i),abs(q)); 


ld=sum(x(:,1:n)’); 
lg=sum(x(:,n+1:w)’); 
z=max(ld,lg); 
Vt-T*z/n; 


ldl=sum(x1(:,1:n)’); 
lgi=sum(x1(:,n+1:w)’); 
21=тах(191,121); 
Vtl=T1*z1/n; 


d=length(find((X-Vt)>0)); 


Pfa=d/S; 
di=length(find((X1-Vt1)>0)); 
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Pfal=d1/S; 


save PFA Pfa Pfal 
quit 
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